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Abstract 

The coming flavour precision era will allow to uncover various patterns of flavour viola- 
tion in different New Physics (NP) scenarios that are presently washed out by hadronic 
and experimental uncertainties. We illustrate this by performing the anatomy of flavour 
violation in the 331 models, based on the gauge group SU(3)c x SU(3)l x U(1)x that 
are among the simplest NP scenarios with new sources of flavour and CP violation. The 
latter originate dominantly through the flavour violating interactions of ordinary quarks 
and leptons with a new heavy Z' gauge boson. After discussing first these models in 
some generality, we present a detailed study of AF = 2 observables and of rare K and 
B meson decays in the specific "/3 = l/\/3" model (to be called 331 model) assuming 
significantly smaller uncertainties in CKM and hadronic parameters than presently avail- 
able. The most prominent roles in our analysis play ek, AM q (q = d,s), the mixing 
induced CP asymmetries S^k s and S^, rare decays B s ^ — > fi + fi~ and in particular 
the CP-asymmetry 3^+ in B s — > [i + ■ As the phenomenology of Z' contributions 
is governed only by Mz> and four new parameters S13, §23, <5i and 62, we identify a 
number of correlations between various observables that differ from those known from 
CMFV models. While, the AF = 2 observables allow still for four oases in the new 
parameter space, we demonstrate how the inclusion of AF = 1 observables can in the 
future identify the optimal oasis for this model. Favouring the inclusive value of \V u b\, 
for 1 TeV < Mz> < 3TeV, the 331 model is in good agreement with the recent data for 
B + — > t + i> t , it removes the ek — S^k s tension present in the SM and the ek — AM Sj ^ 
tension present in CMFV models. Simultaneously, while differing from the SM, it is 
consistent with the present data on B{B S ^ — > and S^. We identify the triple 

correlation B(B S — > — S,^^ — that constitutes an important test of this NP 

scenario. 
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1 Introduction 

The LHCb experiment opened this year a new era of precision flavour physics This 
indirect route to shortest distance scales with the goal to discover new physics (NP) 
is also followed by ATLAS and CMS to be joined later in this decade in particular by 
SuperKEKB, SuperB in Rome and dedicated Kaon experiments NA62 at CERN, KOTO 
at J-Parc and ORKA at Fermilab. 

While, as discussed below and reviewed in [2], there are several anomalies in the 
present data, none of them are conclusive and we should be prepared that NP effects 
in flavour observables will be present at most at the 50% level of the Standard Model 
(SM) predictions and not an order of magnitude as initially expected and hoped for. In 
order to identify NP appearing at this level through flavour violating processes, it will 
be crucial to have both precise experiments and precise theory. 

Indeed in indirect searches for NP through particle-antiparticle mixings in K, B s>< i 
and D meson systems and rare decays of K, B S ( i and D mesons it is crucial to know the 
background to NP: the predictions for various flavour observables within SM. If these 
predictions suffer from large uncertainties also the room left for NP is rather uncertain 
and if a given NP model contains many free parameters, the characteristic flavour violat- 
ing features of this model cannot be transparently seen. They are simply often washed 
out by hadronic and parametric uncertainties even in the presence of accurate data. 
Most prominent examples of this type are the mass differences AM^ S and the parameter 

ex- 
While the sophisticated model independent analyses of NP by professionals like UT- 
fitters [3], CKMfitters [4j, other bayesian |5jj and frequentist analyses [6] and analyses 
using sophisticated Markov-chain Monte Carlos, in particular improved versions of the 
classical Metropolis algorithm [7f[l0| should be appreciated, it is evident from the plots 
presented in these papers that several correlations between various observables are not 
clearly visible when all present hadronic and parametric uncertainties are taken into 
account. 

In our view a complementary and a useful route to the identification of NP is a 
simplified approach in which one investigates how a given NP model would face the 
future more precise data with more precise hadronic and parametric input [^| In this 
manner some of the characteristic features of the particular NP model are not washed 
out and one discovers patterns of flavour violations, in particular correlations between 
different observables, that could distinguish between different NP models. While this 
approach has been already partially followed in the papers reviewed in [2], one recent 
study in this spirit of a concrete model shows that such studies could give interesting 
results. 

Indeed, by proceeding in this manner and analysing the dependence of the allowed 
size of NP contributions in models with U(2) 3 flavour symmetry on the value of \V u b\, 
it was possible to identify the correlation S^k s ~ ~~ Wub\ [14] characteristic for this 
class of models. This correlation was washed out in a much more sophisticated earlier 
analysis in [7] and could not be identified by these authors. Other similar examples in 
the literature could be given here. 

While the more sophisticated analyses listed above are clearly legimitate and possibly 



could appear as the only way to analyze NP models at present, the example in 14 



1 The recent lattice results reviewed in 11 12 and future experimental prospects [T 13 assure us that the 
flavour precision era is just ahead of us. 
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shows that it is also useful to work in a different spirit, even if it is a bit futuristic, and 
investigate how the pattern of flavour violation becomes visible when the uncertainties 
in CKM parameters and hadronic uncertainies are reduced. Such studies are not entirely 
new and can be found occasionally in various proposals, like the one in [13], but here we 
will illustrate them in a concrete model stressing the role of correlations between various 
observables that become visible in this approach. Our strategy is explained in concrete 
terms in Section [H 

One of the important questions in this context is whether the frameworks of con- 
strained Minimal Flavour Violation (CMFV) [15] and the more general framework of 
MFV [16] will be capable of describing the future data. In models of this class, when 
flavour blind phases (FBPs) are absent or set to zero, stringent relations between vari- 



ous observables in K, B® and B® systems are present 17 . Consequently the departures 
from Standard Model (SM) expectations in this class of models in these three meson 
systems are correlated with each other allowing very transparent tests of these simple 
NP scenarios. However, generally these relations can be strongly violated, implying of- 
ten other correlations between observables characteristic for a given NP scenario. Such 
correlations being less sensitive to the model parameters can often allow a transparent 



distinction between various models proposed in the literature [2,18 . In this manner the 
CMFV relations can be considered as standard candles of flavour physics [2]. 

Among the simplest extensions of the SM which go beyond the concept of MFV 
are the so-called 331 models based on the gauge group SU(3)c x SU(S)l x U(1)x 



(331) [19 20 . Here the new sources of flavour and CP violation originate dominantly 
through the flavour violating interactions of ordinary quarks and leptons with a new 
heavy Z' gauge boson. Also one-loop contributions involving new charged gauge bosons 
and new heavy quarks with exotic electric charges can be relevant in certain processes. 

The difference between the 331 models and CMFV models can be transparently seen 
in the formulation of FCNC processes in terms of master one- loop functions [17] . This 
formulation can be used here because the 331 models have the same operator structure 
as the SM and the CMFV models except that the real and universal master functions in 
the latter models become here complex quantities and the property of the universality of 
these functions is lost. Consequently the usual CMFV relations between K, Bj and B s 
systems can be significantly violated. Explicitly, the new functions in the 331 models, 
similarly to the case of the LHT model [21], will be denoted as follows (i = K, d, s): 

Si = \Si\e i9 s, Xi-BlXile&x, Yi = \ Yi \^Y ^ Z^^e^, (1) 

Explicit expressions for these functions in specific 331 models will be given in the course 
of our presentation. 

In 2007 two detailed analyses of flavour violation in the minimal 331 model to be 



defined below have been presented 22 [23]. During the last five years a number of 
changes both in the experiment and theory took place so that a new analysis of this 
model and other 331 models is required. In particular: 



The ek — S^k s tension in the SM has been identified 24-29 . 

The branching ratio for B + — > t + v t has been found to be roughly by a factor of 
two larger than the SM prediction, even though the most recent data seem to be 
consistent with the SM (see Section [5| . 

The new lattice input appears to imply some tension between ek and the mass 
differences AM S)d within the CMFV models |2j. 
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• Most importantly the recent LHCb, ATLAS and CMS data imply SM-like values 
for B{B S —7- and for the mixing induced asymmetry in the B s system, 

whereas in 2007 still large departures from SM predictions for both observables 
were possible. 

The goal of the present paper is to make a fresh look at 331 models and to confront 
their flavour and CP violation structure with recent very improved data on and rare 
decays B s ^ — > fi + [i~ from the LHCb, ATLAS and CMS taking into account the most 
recent lattice input that has a definite impact on various observables, in particular CMFV 
relations. Therefore in addition to and B Sj d — > we address simultaneously the 

£k — S^Ks tension still present in the SM, the mass differences AM S ^ the rare decays 
K + —J- 7r + z/P and Kl — > ir^vv and the tree-level decay B + — > t + v t . We will also consider 
Kl — > B — > X s ^vv and B — > X s j, while leaving the analysis of B — > X s £ + £~ 

and B -)■ K*£ + l~ for the future. 

We should announce already at this stage that in 2012 the minimal 331 model ana- 
lyzed in [2~2~|[23] , the so-called u /3 = \/3" -model, does not imply any longer an interesting 
phenomenology of processes considered by us: new physics (NP) effects are simply very 
small. Therefore, we will dominantly discuss non-minimal 331 models with (3 ^ \^3. 
While a significant part of our presentation will apply to arbitrary (3, we will present a 
detailed phenomenology of flavour observables in a specific 331 model with f3 = 1/^/3 
to be named 331 model in what follows. This model, while being rather economical, 
implies much more interesting phenomenology than the minimal model. In particular, 
while consistent with the present data, among rare decays, the largest deviations from 
SM expectations are found in B s ^ — > /i + branching ratios and the CP asymmetries 
S s t _mB sd -> , recently introduced in 



30 
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Our paper is organized as follows. In Section 2] we recall in some details basic fea- 
tures of 331 models that are relevant for understanding of our analysis. In particular, we 
derive all Feynman rules for this class of models that we need for our phenomenology. 
In Section [3] we present the formulae for the effective Hamiltonians governing particle- 
antiparticle mixings K° — K° and Bd jS — Bd,s that in addition to the usual SM box 
diagrams receive tree-level contributions from Z' exchanges and new box diagram con- 
tributions. We also give a compendium of formulae relevant for numerical analysis of 
AF = 2 observables. In Section [4] the effective Hamiltonians for s — > dvu, b —> qvv, 
b — > ql + £~ (q = d, s) and b — > sj transitions are given. In particular we present an im- 
proved QCD analysis of Z' contribution to B — > X s j, demonstrating that it is negligibly 
small compared to the SM contribution. In Section [5] we calculate the most interest- 
ing rare decay branching ratios in the K and B meson systems, including those for the 
processes K + — > tt + vv, Kl — > tt°vv, B — > X s ^vv, B s c i — > [J. + n~, B + — > t + v t and 
Kl — > / i+ A i • In Section [6] we recall the most important CMFV relations with the goal 
to find out how these relations are violated in the 331 models. In Section[7]we emphasize 
stringent relations between various observables in the 331 models that originate in the 
dominance of NP contributions by the tree- level Z' exchanges. In fact even without a 
detailed numerical analysis, it will be possible to deduce from these relations the pattern 
of deviations from the SM expectations predicted by the 331 model. In Section [8] we 
present first our strategy for numerical analysis and summarize the recent LHCb data 
and new lattice input. Subsequently, we demonstrate how already AF = 2 observables 
allow to identify four allowed oases in the new parameter space. We demonstrate how 
the inclusion of AF = 1 observables can in the future select the optimal oasis for this 
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model. In this context the branching ratios B{B Sj( i — > ), the CP-asymmetry S^ +fi _ 
in B s — > fi + fi~ and S^^, when considered simultaneously, play the prominent roles. The 
correlations between these observables in the 331 model should allow to monitor how 
this model will face the data in the coming years. We also analyze how these correlations 
differ from the ones present in the CMFV models. The highlights of our analysis are 
listed in Section [9j There we also make a comparison of our findings with those known 
from other simple models. Few technicalities are collected in Appendices. 



2 Flavour Structure of the 331 Models 
2.1 Preliminaries 

Among the many scenarios that have been proposed to extend the Standard Model (SM), 
one with rather appealing features is that based on the gauge group SU(3)c x SU(3)l x 



U(l)x, originally developed in 19,20 . This group is spontaneously broken to the SM 
group SU(2)l x U(l)x, subsequently broken down to U{1)q. Therefore, the 331 model 
has an extended Higgs sector [32] , with the first symmetry breaking occurring at a scale 
much larger than the electroweak scale. 

One of the nice features of this model is that the requirement of anomaly cancelation, 
together with that of asymptotic freedom of QCD, constrains the number of generations 
to be equal to the number of colours, thus providing an explanation for the existence of 
3 generations, long sought for in the SM. Requirement of anomaly cancelation has also 
consequences on the transformation properties of fermions. In fact, as a first consequence 
of the extension of the SM gauge group SU(2)l to SU(3)l, one has that left-handed 
fermions transform as triplets (or antitriplets) under the action of SU(3)l. In order to 
have an anomaly-free theory, the number of triplets should be equal to the number of 
antitriplets. A possible choice is to assume that the three lepton generations transform 
as antitriplets, so that, taking into account the three colour possibilities for the quarks, 
the number of quark triplets should be equal to the number of antitriplets plus one. 
Hence two quark generations should transform as triplets, one as an antitriplet. The 
common choice of having the third generation with different transformation properties 
might be at the origin of the large mass of the top quark with respect to the other 
quarks. However, in contrast to the SM, where anomaly is canceled for each generation 
of fermions, in the 331 model the cancelation is fulfilled only when all the generations 
are considered. 

As for the content of the (anti)triplets, the upper two components coincide with 
the known SM fermions. In the case of quarks, the third members of the generations 
correspond to new heavy quarks, while for leptons different choices are possible. The 
minimal version of the model makes the choice that no new leptons are introduced, so 
that the third component of the lepton antitriplet is chosen as the conjugate of the 
charged lepton (in practice, all the three possible helicity states for the SM leptons are 

contained). In this minimal version, among the three new heavy quarks two have charge 

4 5. 
— — and one — . Other variants are possible. In particular, in the lepton sector, another 

possibility could be that of having new heavy neutrinos as third components of the 
lepton antitriplets |33| . Indeed in the present analysis we choose a variant of this latter 
kind that has phenomenological implications that are more interesting than found in the 
minimal model. 
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The extension of the SM gauge group SU(2)l to SU(3)l also implies the existence 
of 5 new gauge bosons. These are a new neutral boson (Z') plus other four that might 
be charged depending on the variant of the model that one chooses. We denote them 
generically as V^ v and Y^® Y . As we shall see when describing the model in details, the 
new bosons V^ v and Y^ Y couple two SM leptons, while this is not possible for SM 
quarks, which can only be coupled to new heavy quarks by means of these gauge bosons. 
Therefore the latter are assigned lepton number L = =p2, but carry no lepton generation 
number so that the lepton generation number can be violated and indeed processes like 
(j, — > e7 or n — > 3e might proceed due to such new gauge boson mediation. Therefore, 
bounds on this model could come from the MEG experiment [34] the physics programme 
of which is devoted to the search for the lepton flavour violating process \i — > ej. Detailed 



studies of the lepton sector of the model have been performed in [35 36 , in particular 
in [35] a mechanism to suppress large CP violation in this sector as well as large leptonic 
electric dipole moments was found. 

The new neutral gauge boson Z' also introduces new features with respect to the 
SM as well as with respect to other NP scenarios that also contain a new neutral Z', as 
for example grand unified theories (GUTs). In fact, in contrast to most of them, the Z' 
mass in some 331 models cannot be arbitrary, but is limited from above. R This occurs 
because the U(l)x gauge coupling gx and the SU(3)l one g are related by: 

g\ _ 6sin 2 w 



g 2 l-(l + /3 2 )sin 2 w 



(2) 



where 9\y is Weinberg angle and (3 is a parameter of the model (see next subsection) 

that defines its variants. Therefore, considering for example the minimal model with 

1 

/3 = \/3, one should require sin 2 9w(Mz>) < -■ Since we know that sin 2 9\v(Mz) — 0.23, 

the evolution up to Mz' constrains Mz> below a few TeV [38] . On the other hand for 
(3 = l/\/3, as dominatly considered in our paper, this bound does not apply. 

Also lower bounds on Mz> were discussed in literature, stemming from study of [i 
decay within the model 1 35, 38 ,39], from the analysis of S, T, U parameters [33] and from 
comparison with experimental results from CDF Collaboration for possible Z' decays to 
e + e~ or fi + ji~ |40| . In principle those bounds are model dependent, since they depend 
also on the entries of the unitarity matrices that transform the quark gauge eigenstates 
into mass eigenstates. However, taking into account the results of the aforementioned 
analyses, we adopt as working range 1 TeV < Mz> < 3 TeV, which also shows that 
present experiments, mainly those at the LHC, have the possibility to test the model in 
the nearby future. Along the same lines, a possible mixing angle between the SM Z and 



the new Z' has been considered, with the result that Qz-Z' < 0(10 ) 40 . In what 
follows we will neglect this mixing. 

The main novelty with respect to the SM related to the Z' is the existence of tree level 
flavour changing neutral currents (FCNC). However, these arise only in the quark sector, 
since the universality of the coupling of the Z' to leptons guarantees that no FCNC show 
up in this case. Furthermore, universality is also realized in the Z' couplings to right- 
handed quarks, so that the new FCNC are purely left-handed. The relevant couplings 
in this case are anyway small enough to make such tree level transition not larger than 
the corresponding loop induced SM contribution. Many processes induced by FCNC 



2 Another similar case is for example represented by heterotic string models with supergravity mediated 



supersymmetry breaking 37 where also the Z' mass cannot be larger than around a TeV. 
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have been considered in the framework of the 331 model. Mixing of neutral mesons as 
well as a number of rare K and B^ s decays have been considered in 39 , 4lJ . Some of 



the results have been updated in 22 , with particular attention to the neutral meson 
mixing and the study of possible new CP violating effects. The main result contained 
there is that the most sensitive quantity is the phase of B s mixing /3 S , as determined from 
B s — > J/t/j(f>. This is partly a consequence of the expectation that purely hadronic decays 
are potentially more sensitive to the Z' contributions since in addition to being flavour 
violating, the flavour conserving couplings of Z' to quarks is enhanced with respect to its 
coupling to leptons. We shall see this explicitly below. Finally, b — > sj was considered 
in [42] and reanalysed in |23| . Notice that this mode is loop-induced also in the 331 
model. 

As a final remark, we mention that the 331 model naturally possesses the Peccei- 
Quinn symmetry [43], solving the strong CP problem [44]. 



2.2 The Model 
2.2.1 Gauge Bosons 

In this section we provide a description of the 331 model in general terms and some of 
its variants parametrized by j3. 

Leaving aside the case of SU(3)c, that presents no differences with respect to the 

A a 

SM, the generators of the gauge group are: T a = — for SU(3)l (A a being the Gell- 

Mann matrices and a = 1, ... 8) plus T 9 which generates U(l)x- T 9 can be conveniently 
1 

defined as T y = —=, the unit matrix being 3x3, when acting on triplets and just 1 
V6 

when acting on singlets. With this definition, the generators satisfy the same relation as 
the eight generators of SU(3): Tr[T a T b ] = ^ . In the case of antitriplets, transforming 

as 3, the generators are T a = — (T a ) T , T meaning transpose. 

The gauge bosons of the group SU{3)l are denoted as W a , a = 1, ... 8 and transform 
according to the adjoint representation of the group, while the gauge boson of U(l)x is 
denoted by X. As a consequence, the covariant derivatives acting on the various fields 
read as follows: 

• triplet ipL: D^ L = Of/tpL -igW*T a ij) L - ig x X X„T 9 ^ L ; 

• antitriplet 4> L : D^ L = d^ L + i g W*{T a ) T $ L - ig x X X^T 9 ^ L ; 

• singlet tp R : D^ R = d^ R - ig x X A^T^R- 

In the previous formulae, g, gx denote coupling constants and X is the new quantum 
number corresponding to U{l)x- 

The matrix W a T a can be rearranged as follows: 



W» = W«T a = 1 



( W 3 + _^ W 8 ^yQr \ 

V2W- -W* + j= 3 W* V2V? V 
V2Y^ Qy V2V,r Qv -^W»y 



(3) 
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where the following definitions have been introduced: 

Y**y = ^{WiTiWl), (5) 

Vf Qv = ^{W^iWl). (6) 

W 3 and W s correspond to neutral gauge bosons; to the charged gauge bosons of the 
SM, while V^ v and Y^ Y might be charged or neutral. In order to assess this precisely, 
we have to define the electric charge matrix. In analogy with the SM, we introduce the 
following operator: 

Qw = f 3 + | (7) 



Y 
~2 



pf 8 + Xt . (8) 



The constant /3 that has been just introduced plays a key role in the model: according 
to its value various versions of the model have been proposed, which significantly differ 
from each other. In particular, only for some values of (3 the gauge bosons turn out 
to have integer charges. In the case of the minimal 331 model, the choice j3 = v3 is 
adopted. 

Even though we shall recall the features of the minimal model, we shall instead 

consider the variant with f3 = which is phenomenologically more interesting than 

v3 

the (3 = V3 case. Moreover, in contrast to the minimal model, it does not introduce 
exotic electric charges of fermions in order to cancel the anomalies. We also note that 
|/3 1 = 1/V3 is the smallest value of this parameter for which the heavy gauge bosons 
have integer electric charges. 

In (7j8), the hat reminds us that we are dealing with operators and their action 



should be distinguished from the simple multiplication with the corresponding matrices. 
Let us consider in turn the various cases: 

• Triplets: Qw = Qw Plet ipL with 



^triplet 



(\ + ^ + x 



\ + in + x I (<):i 

~73 +X 



Antitriplets: Qw 4>L = Q^ tltr%plet il)L with 

Qantitriplet 
W ~ 

\ 



• Singlets Q w 4>R = Q^ 9et ipR with Q 3 ^ 9 ^ = X. 




(10) 
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Gauge bosons: QwW^ = [Q^-, = Q^pW^ so that the charges of the entries 



in (|3j) can be read from the matrix Q^p (GB standing for gauge bosons): 



n GB _ 



1j -p4 \ (ii) 



2 



We can see that there are three neutral gauge bosons, two singly charged ones, that 

can be identified with the SM bosons, and other four bosons. Choosing /3 = — -= 

V3 

these are two singly charged bosons and two additional neutral bosons V°, V°. For 
comparison, in the minimal model one finds instead two new singly charged bosons 
and two new doubly charged ones 

After these general discussion we now consider the particle content of the model, 
that is fermions, Higgs particles and gauge bosons. The main target in this step is 
the derivation of the interaction Lagrangian in the mass eigenstate basis for all fields 
involved. This will allow to find directly the relevant Feynman rules. We begin with 
fermions. 






2.2.2 Fermions 

As already stated, left-handed fermions fit in (anti)triplets, while right-handed ones are 
singlets under the action of SU(3)l- In the case of the quarks, the first two generations 
fill the two upper components of a triplet, while the third one fills those of an antitriplet; 
the third member of the quark (anti)triplet is a new heavy fermion. Therefore we have: 



(12) 



The corresponding right handed quarks are singlets. Notice that in the third generation of 
quarks a different order is chosen with respect to the SM. This stems from the requirement 
of the anomaly cancelation, as discussed at the beginning of this section. The sign of 
the fields in the antitriplet is fixed requiring that the SM couplings are reproduced. 

The case of leptons is more involved. As already stated, left-handed leptons fit in 
antitriplets, and again the SM fields fill the upper two components of the antitriplets. 
The third member of the antitriplet is chosen in a different manner depending on the 
considered variant of the model. In the minimal version it is just the conjugate of the 
charged lepton in the same antitriplet, while in a generic version of the model it could be 
a new heavy lepton. In the latter case, one needs to introduce the right handed charged 
leptons as singlets. Therefore, the antitriplets look as follows. 





(13) 



with E e = e c , En = /u c , E T = t c only in the minimal version. Hence, in the minimal 
model one can avoid the introduction of new leptons as well as right-handed leptons 
in separate multiplets. Right-handed neutrinos are not compatible with the minimal 
model. In the case of /3 = I/a/3 we have E e = v%, E^ = E T = 
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arbitrary 




fermions 


Qw 


X 


Qw 


X 




2 
3 


1 p 

6 2^ 


2 
3 





UR, Cr 


Z 
3 


z 

i s 


z 

3 


z 

3 


dL, s L 


1 

~~ 3 


6 2^ 


1 

~~ 3 





j 

UR, Sr 


1 
1 

~~ 3 


1 

_ 3 


1 
1 

~~ 3 


1 

~~ 3 


D l , S l 


1 V3P 
6 2 


1 /8 
6 2V3 


1 

3 





Dr, Sr 


i Vsp 

6 2 


1 \/5j8 

6 2 


1 

3 


1 

3 


b L 


1 
3 


6 ^ 2^ 


1 

3 


1 

3 


t L 


2 
3 


1 _L _£_ 

6^2^ 


2 
3 


1 

3 


b R 


1 

~3 


1 

J 


1 

— 


1 

—-3- 


t R 


z 

3 


Z 
3 


Z 

3 


z 
3 


T L 


i _|_ Vzp 

6 ' 2 


I -1 — £_ 


2 
3 


1 

3 


Tr 


1 , V3/3 
6 ' 2 


1 , v^/? 
6 2 


2 
3 


2 
3 


h 


-1 


— - + -2- 

2^2^ 


-1 


1 

3 


t R 


-1 


-1 


-1 


-1 


VlL 





2^2^ 





1 

3 


VlR 














EgL 


1 , V3/3 

2 ' 2 


2^2^ 





1 

3 


Em 


1 , Vsp 

2 ' 2 


1 , v^/? 

2 ' 2 









Table 1: Electric charges and X quantum numbers of the fermions in the 331 model for 
arbitrary (3 and for /3 — 4=. In the case of leptons I = e, /i, r. 



n-n i j. ■ i r n r i_ j r ^triplet ^antitriplet 1 ^sinqlet 

ihe electric charge of all termions can be read trom Q w , Q w and Q w . 

Furthermore, imposing that the SM fermions have the usual charge, the value of X for 
such fields can be fixed. Due to gauge invariance, X is the same for all the members of a 
given (anti)triplet. Therefore, fixing for example the charge of the ul quark 2+^3+-^ = 

| allows to assign the value X(u£) = \ — which in turn fixes the charge of the new 

heavy quark Dl in the same triplet to = — + X = g — ^y^, while correctly 
reproducing the dL quark charge. Finally, for right-handed fields one has Qw = X, 
Qw being of course the same as for left-handed fields. In Table [T] we collect the electric 
charges of all the fermions and their X quantum numbers, both for arbitrary (5 and for 

8 = — =. In the latter case we see that no exotic charges appear. 
2.2.3 Higgs sector 

Symmetry breaking is accomplished in two steps. At first one has to recover the SM as a 
low energy effective theory deriving from the 331 model after the spontaneous symmetry 
breaking SU(3)l X U(l)x — > SU(2)l X U(l)x at a given scale, that has to be properly 
identified with a vev of a Higgs field. Denoting with < <p\ > such a vev, in this step it 



2 Flavour Structure of the 331 Models 



10 



should be required that 

f 1 < fa >= f 2 < fa >= f 3 < fa >= (fit 8 + XI) < fa >= , (14) 

while all the other generators, when acting on the vacuum, give a non vanishing result. 

In the second step, the familiar SM symmetry breaking SU(2)i x U(l)x U(1)q 
should be realized at the electroweak (EW) scale, the corresponding Higgs vev being 
denoted by < fa >■ Now the only generator that leaves the vacuum invariant should be 
Qw- 

Qw < fa >= . (15) 

Since the final goal is to give masses to the fermions and all the gauge bosons except 
the photon, the procedure to be followed is still the building of a Lagrangian for the 
Higgs fields containing the coupling to gauge bosons through the covariant derivative 
and a Yukawa Lagrangian coupling such fields to fermions, in such a manner that the 
symmetries of the model are preserved. In particular, the X conservation should be 
guaranteed. It turns out that invariant terms involving a fermion triplet and a singlet 
can be built introducing a Higgs field that transforms either as a triplet or as a sextet. 



However, in the case of the sextet, requiring that the conditions (14) are satisfied con- 
strains the parameter f3 to specific values (different for quarks and leptons) that are not 
compatible with the scenario considered in this paper. Therefore, we are left with the 
choice of a Higgs triplet, that we denote as x- I n this case, the first three conditions in 



(14) constrain the vev of x to the following structure: 

< X >=4= I | • (16) 



quantum number X for \ to the value X x 




Furthermore, imposing the last condition in (14), i.e. (/3T 8 + Xt) < x >= 0, fixes the 

T 

We now turn to the second step in the spontaneous symmetry breaking. Again, the 
Higgs field can transform under the action of SU(3)l either as a triplet or as a sextet. 
Let us start with the first possibility. Requiring that all the generators except Qw are 
broken implies that one can choose two possible Higgs triplets, with the following vevs: 



<p>=— [ v \ , <^>= o . (17) 





The last condition: Qw < p >= Qw < V >= finally fixes: 

xp= \~w* Xn= ~\~w* (18) 

The case of the sextet can be carried out following the same steps and recalling 
that the action of the SU(3) matrices T a on a generic sextet S is defined as: T a S = 
ST a + (T a ) T S. The result is that there are four possible independent sextets. A usual 
choice, that we adopt in this paper as well, is to keep only the following Higgs sextet: 

1 f° °\ 1 B 

'--A::: * x ^ + ^- (19) 
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The X value reported above has been fixed as done in the case of the Higgs triplets. 

Notice that, since the EW symmetry breaking should occur at a scale much lower 
than the scale at which SU (3) l is broken down to SU(2)l, the following hierarchy among 
the vevs should hold: «>», v' , w. 

The charge of the fields in the three Higgs triplets X: Pi V can he read from Q tri P let 
by inserting the corresponding value of X. As for the Higgs sextet S, the electric charge 
matrix is: 



Q se*tet = \ A + 2X _ X + A + 2X _1 _ J» + 2X 



± .. •> v J_ _l •> v I - _L T _ v v 
^=-r^ -^t^t^a -j-j^+2X . (20) 

It turns out that the electric charges for complete set of Higgs fields can be expressed 
in terms of the following /3-dependent quantities Qa = \ + Qb = — \ + and 

Qc = -§ + ^: 

Qa 0^ 

Q x = | Q B | , Q p = | | , Q v = \ -1 J . (21) 


Qa Qb 

Qs= | QbQc -1 | • (22) 
-1 

For /3 = -^= one has: = 1, <5_B = and Qc = —1 (for comparison, in the minimal 
v3 

model one finds Qa = 2, Q# = 1 and Qc = 0). Since the sextet does not enter in 
the quark-Higgs interactions, we shall only concentrate on the three triplets and, on the 
basis of their electric charges, we adopt the following notation: 









X = \XB , p = P u , V = tT ■ (23) 



The neutral fields x°, p° and if can be decomposed into their real and imaginary 
parts: 

X° = 4(& + iCx). (24) 



V2 

P° = ^(Cp + *Cp) , <£p>=^, (25) 
if = -L(^ + iC„) , <£„>=«', (26) 

with the ratios of the vevs defining three angles that govern the mixing among the fields 
in this rich Higgs sector: 

v' 2 v 2 v' 2 

sin 2 f} v , v = -;r- — H , sin 2 VU = — - , sin 2 [} v , u = — . (27) 

yZ _|_ yli yZ _|_ yZ yZ _|_ ylZ 
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Introducing the Higgs potential: 

Vmggsip, V, X) = l4{p ] p) + Ai(pV) 2 + plirfv) + A 2 (r/ t r/) 2 + Mx) + Hx j xf 

+ Ai 2 (p t p)(?? t »7) + Ai 3 (p t p)(x t x) + ^23(v t v){x t x) + 

+ Ai 2 (p t r ? )(r ? t p) + Ai 3 (p t x)(x t p) + M^xXx^) + 

+ V2f 1 (e ijk p i r r i x k + h.c.) (28) 

and minimizing it with respect to u, v, v' one can express [if, i = 1, 2, 3 in terms of the 
other parameters in VHig gs (p,T]iX)- 

Higgs fields mix, giving rise to the mass eigenstates listed below. 

• Singly charged states having \L = 0|: 

^ "\ _ ( COS ~ sin Pv'v \ ( P ± 



Hq J \sm(3 v > v cosP v f v J \v ± J 
with M,± = and M 2 + = . R a — , where we have exploited the hierarchy 

<P W H£ sin P v , v cos /?„,„ ' ^ J 

u v, v'. 

States with L = 2 and charge Qb: 



<l>v\ = ( cos ^3 -sin/3 13 \ / p B 
with M,s = and M 2 B = /i^tan + u 2 — — and /3i 3 = f3 v 

V V n 2 2 2 

States with charge Qa- 

I) y \ — ( cos /^23 - sin /3 23 \ ( rj A 



(30) 



(31) 





V sin /523 cos/3 23 )\x A 

with M^a = and = fiucot/3 v > v + u 2 ^- and /3 13 = ^ - 

• Neutral pseudoscalar states: 

cos d sin /Sy'u 

sin/3„/„ cos/Vu | | ^ x | (32) 
0-1 

with M? = M? =0 and M? = — . The real parts of the neutral 

412 2sin/V*,cos/V„ P 

Higgs fields gives raise to scalar states, with a complicated mixing pattern. Quarks and 
gauge bosons couple to the various Higgs fields defined by the above mixing. 

2.2.4 Masses of the gauge bosons. 

Having fixed the quantum numbers of the various Higgs fields, we can proceed to derive 
the masses of the gauge bosons. These stem from the covariant derivative in the Higgs 
Lagrangian: 

L mggs = {D„ X )\D» X ) + (IV) f (£"V) + (D^)\D»n) + Tr[(D,S)\D»S)] . (33) 
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The masses of the W ± , V ±Qv , Y ±Qy bosons read: 



M w± 



M, 



9 2 2 
9 2 2 

1 U 



M, 



Y±Qy 



21 



(34) 
(35) 
(36) 

v 2 — w 2 . The remaining 
neutral gauge bosons mix with each other. In a first step, the two neutral gauge bosons 
and mix, giving rise to the two bosons and Z 1 ^. The mixing angle is denoted 



9 2 2 
4 U 



v 2 _ v\ 
2^ + 2^2 



v_ v+ 
2u 2 + 2u 2 



where we have defined: v\ 



v + v + w and v_ 



/2 



by 6I331 and is given by sin 6*331 



9 



and mix with W 3 . However, 



taking into account the hierarchy u 3> v, v' , w, it turns out that the mixing involves 
only Bfj, and W 3 finally giving a massless photon and the to be identified with 
the SM Z. Introducing a new constant 

1 6 



9\ 



1 +" g 2 



(37) 



it turns out that the mixing angle that defines the rotation matrix providing the mass 
eigenstates A, Z starting from W 3 , B is given by 



sm9\Y = 



9Y 



^9 2 +9y 



(38) 



Therefore, in this model the following relation between the two mixing angles 6*331 and 
6w holds: 

cos 6*331 = (3 tan $w ■ 



Thus the mixing of these neutral gauge bosons can be summarized as follows: 



The resulting masses for the neutral gauge bosons are: 



Ma 



M 2 Z 



M 2 Z , 



, 



9 



4cos6V + ' 

g 2 u 2 cos 2 0w 
3[1- (l + (3 2 ) sin 2 6 W ] 



1 + 



4u 2 



1 + 3/3 2 



2 sin 4 Ow \ \/3/3 v 2 sin 2 6 W 



cos 4 0y/ 



+ 



2 u 2 cos 2 6w 



(39) 
(40) 

(41) 
(42) 

(43) 



Considering again that u 3> v, v' , w one is lead to the final expressions for all the 
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gauge boson masses: 

M W± = 9 -T v li ( 1 1 ) 



4 

a 2 

r2 _ A/f2 _ i» „,2 



M' v±Qv = My ±Qy = »-u z , (45) 



4 

r2 



Mi = , (46) 
4 cos 6% 



r2 _ g V cos 2 6> w 

3[1- (l + /3 2 )sin 2 W ] 



M Z> = oM „ , 2 „ n • ( 4 8) 



M 2 

Therefore, as in the SM, one has that = cos 6w, which allows to identihy the angle 
Ow with the Weinberg angle. 

2.2.5 Yukawa Lagrangian and interaction Lagrangian terms in their 
final form. 

The next step is the definition of the quark mass eigenstates. The pattern is similar to 
the 2HDM or SUSY one, with up-type quarks getting mass from the coupling to the 
Higgs triplet p and down- type quarks from the coupling to the Higgs triplet r\. The new 
quarks D, S, T get mass from the coupling to the Higgs triplet X- The most general 
Yukawa Lagrangian coupling left and right handed quarks to Higgs fields in a gauge 
invariant way is 

Lyuk = ^i, a QiPda,R + ^aQzrf d a ,R 
+ KaQiVUa,R + A^ o g 3 P*«a,fl 

+ *ijQiXJj,R + H,zQzX*T R + h.c. (49) 

where Qi, i = 1,2, are the left-handed triplets of the first two quark generations while 
Qz is the corresponding one of the third generation; = D, S; a = 1, 2, 3 with 

d (i,2,3)R = d R, s R , b R and u^ 3 )r = u R , c R , t R . 

Upon rotation through unitary matrices one defines the mass eigenstates: 





A 




(50) 



where the rotation matrices are unitary 

U\U L = U L U[ = VIV L = V L VI = 1 (51) 
and their matrix elements will be denoted as follows 

Vij = (V L )ij, Uij = (U L )ij. (52) 

They will be crucial ingredients in the interaction Lagrangians below. In the following 
we shall omit the prime in the notation for these fields, understanding that we are 
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considering the mass eigenstates. It can be noticed that, among the new quarks, T 
cannot mix, since its charge is different from that of the other new quarks. 

As for the D , S quarks, in principle they could also mix; however, as discussed 
in 22 , 39 , a convenient choice is that of taking them simultaneously as gauge and mass 
eigenstates. Indeed the mass matrix Mj of D, S, T is block-diagonal, since T cannot 
mix. With two block-diagonal rotation matrices Wl and Wr, Mj can be diagonalized. 
Then Wr can be removed by a redefinition of right-handed heavy quark fields and Wl 



can be absorbed into the definition of Ul and Vl- Below in Eq. (64) we will choose a 
parametrization for Vl that is compatible with this redefinition as discussed in [22] . 

As in the SM, the rotation of the quark fields affects the charged currents. In complete 
analogy with the SM, it is possible to define a CKM matrix as 

Vckm = Ulv L , (53) 

which weights the transitions between up-type and down-type SM quarks belonging to 
different generations mediated by the W^. 

However, in the present model also neutral currents mediated by the Z' are affected by 
the quark mixing. This happens because the couplings of the Z' to the third generation 
of quarks differ from the ones to the first two generations resulting in tree-level FCNC 
mediated by Z' . In order to see this explicitly, we look at Z' couplings to SM quarks 

C z ' = J^Z'v , (54) 
Jfi = ul7v U l a U L u L + dLlnVl I a I V L d L , (55) 



where the couplings a and b depend on the chosen (3. Separating the flavour conserving 
and flavour changing parts we find 



+ u Llll ul U L u L + d Llfl v[ 





(56) 



which shows that the different treatment of the third generation in 331 models generates 
FCNCs at tree level. 

We are now ready to list the interaction Lagrangian densities (multiplied by i, in 
order to directly have the Feynman rules) that give the couplings of fermions to various 
gauge bosons. We list first the terms containing the currents generalizing the analogous 
SM ones and then those containing the currents mediated by the new gauge bosons. 



• charged current mediated by 



t=e,fi,T 



+ i-%= I (u L c L t L ) I^Vckm I s L I W + » + h.c. I (58) 



V2 
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• e.m. current 

iL] nt = J2iQeh»fA», (59) 
/ 

/ is a generic lepton or quark field of charge Q in units of the positron charge 
e = g sin Ow- 

• neutral current mediated by Z boson 



2cw 

X ]ytLl^tL - (1 - 2s^fcA + 1s 2 w Ir1^ r + (1 - V3f3)slrE a ^E e 



q u =u,c,t ^ ' 
+ E [ ( _1 + \ S W ) VdLl^qdL + ^S^qdRl^dR 



+ 2 (-^ + ^) ^ + ^] -2^ + ^) S ^ T } ( 6 °) 

where we adopted the notation sw = sinOw, cw = cos Ow- As in the SM, for 
ordinary quarks this current contains a piece proportional to the e.m. current plus 
a term in which only left-handed fermions are involved. In the case of the new 
quarks, only the contribution proportional to the e.m. current is present. 

current mediated by V ±< ^ v 

9 



Lf nt = ij= X [viL^E eL V- Q v» + E eLlll v eL V +Q ^ 

■ 9 

'T2 



£=e,fi,r 

+ E [{^MQdl^V^ + {QDL)MqdL)iVijV-Qv*\ 



=1,2,3; j=l,2 



+ X [ f LlMuL)iUziV + Qv» + (q^a^TLV^V-^] } (61) 



=1,2,3 



where q u i i = 1,2,3 correspond to u,c,t; qdi i = 1,2,3 correspond to d,s,b; Qr>% 
% = 1,2 correspond to D, S. The elements of the matrices Ul and Vl have been 
denoted by Uij and Vij respectively. 

current mediated by Y ±< ^ Y 

iL Y nt = -ij= X [EiL7»l L Y + ^ + £ L ^E eL Y-Q^] 

+ E [(quL\l^QDL) U* 3 Y+QY ^ + (QDL^iquL^UijY-^] 

V 1 i=l,2,3;j=l,2 

- E [TLlMdL^vuY^^ + iqdLhl^Lvl^^]} (62) 
i=l,2,3 
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• neutral current mediated by Z' 

-z> , gz" 1 



2v / 3c l y^l-(l + /3 2 )s^ 

E {[i-(i+V3^) 

k £=e,n,T 

- [2 - (2 - \/3>(l - >/3/3)) s^] EurfpEtL - V3(l - y/SpWs^EuapEui} 

ij=l,2,3 V3 

+ [-! + (! + (MdL)iln{qdL)j$ij + ^{qdL)afi{qdL)jvliV Zj 

4 2 1 



r r -2\/3 + /3(l - 3^/3/3) 2 n , = , ^ , 

+ E {[ 2 + — ^ —4] {QDLia^iQDL) 



1=1,2 



2 v / 3 + /3(l-3 v / 3/3) 

g(l - 3^3/3) 2 = 
H s w (Q DR )a^(Q DR ) 



r 2^ + ^(1 + 3^/3) 2 1 - 5(1 + 3^,- t^IU^ 
+ L ~ 2 + -y| s w \T Ll ^T L + ^= s w T R j^T R j U63) 

It can be noticed that while the current mediated by the SM bosons depends on 
the CKM matrix elements, the currents mediated by the new gauge bosons V, Y and Z' 
depend on the elements of the two rotation matrices Ul and Vl introduced above. Due 
to the relation U\Vl = Vckm one can choose Vl and Vqkm as independent matrices. 
In the following we adopt the parametrization of the matrix Vl reported below [22 [^J 

/ C12C13 Sl 2 C 23 e^ 3 - C 12 Sl3S23e l( - Sl ~ S ' 2) Ci2C23Sl3e lSl + S!2S23e l( - S2+53) 

V L = -ci 3 si 2 e" l<53 C12C23 + 5i2Si3S23e j(<5l_52_53) -si2Si3C23e i(l5l ~ <53) - ci 2 s 23 e l<52 

\ -S\3e~ lSl -Cl3S23^ lS ' 2 C13C23 

(64) 

One advantage of this parametrization is that the elements of the third line which will 
appear in the Z' contributions to FCNC observables have a small number of parame- 



ters, e.g. si 2 , C12 and the CP violating phase £3 do not occur. As pointed out in 22 
this parametrization is compatible with the choice to treat D and S both as mass and 
interaction eigentstates. 



A similar parametrization had been proposed in [39 with the difference that the weak phases were ne- 



glected. 
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3 AF = 2 Transitions 



3.1 Preliminaries 

In order to have transparent formulae for various observables we introduce a compact 
notation for the couplings presented above. First of all we denote the Z' couplings to 
down-quarks as follows: 



iC L {Z') = i Af(Z')(s^P L d) + AZWfrfPLd) + A b L s (Z')(b^P LS 



Z' 



(65) 



with the first upper index denoting outgoing quark and the second incoming one. Con- 
sequently 

A^(Z') = {Ai{Z'))\ (66) 

The expressions for A^(Z') in terms of the fundamental parameters of the 331 models 
considered are collected in Appendix [A] There the Z' couplings to leptons, like A^(Z') 
and A^ /J R (Z / ), are defined in analogy to (65). 



3.2 Standard Model Results 

We begin the presentation of AF = 2 transitions by summarizing the known results 
within the SM. SM contributions to the off-diagonal elements M[ 2 in the neutral K and 
B q meson mass matrices are given as follows 



M- 



K\* 
12;SM 



Gp 



12tt 2 

( M i 2 )sm = ^2 F kBB d m Bd M 



FkBkttikMw [\ 2 c r}iSo{x c ) + XtV2S (x t ) + 2X c X t r] 3 So( 

) ,(67) 



AS 9) ) VB S (x t ) 



where X{ = mf/Myy and 



A 



(K) 



V£ s V id , 



A 



(9) 



VttV tg 



(68) 



(69) 



with Vij being the elements of the CKM matrix. Here, So(xi) and Sq(x c , x t ) are one-loop 
box functions for which explicit expressions are given e. g. in [45]. The factors r/j are QCD 



51 



52 



corrections evaluated at the NLO level in |46 - 50 . For r/i and 773 also NNLO corrections 
have been calculated 



Finally Bk and B B are the well-known non-perturbative 



factors. 

In the SM only a single operator 

QX^(K) = {s llx P L d) (jrfP L <t) 



QY LL (B q ) 



{h^P L q) (bj»P L q) (70) 



contributes to and Mf 2 {q = d,s), respectively. Moreover flavour and CP violation 
is governed totally by the CKM matrix. 

In the 331 model the operator structure remains unchanged, which does not increase 
the hadronic uncertainties relatively to the SM. But there are new flavour violating and 
CP-violating interactions originating dominantly in tree-level Z' exchanges shown in 
Fig. [TJ Also box diagrams with new charged and neutral gauge bosons and new quark 
exchanges can contribute and we comment on their importance soon. 
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Figure 1: Tree-level flavour changing Z' contribution to K° — K° mixing (the diagram rotated 
with 90° also exists). 



The effect of these new contributions can be summarized by replacing the flavour 
independent So(xt) in the SM formulae by the functions Si (i = K,d,s): 



Si = So(x t ) + ASi(Z') + AS; (Box) = \Si\e ie s . 



(71) 



The important new property is the flavour dependence in these functions and the fact 
that they carry new complex phases. In what follows we will give the formulae for the 



new contributions in ( 71 ) . 



3.3 Tree Level Z r Contributions 

We begin our discussion with the tree level Z' exchanges contributing to AF = 2 tran- 
sitions in Fig. [TJ Defining 

2 a Gf 



a 



we find 



9SM 



^2 2vr sin 2 W 



(72) 



A5,(Z') 



A b «(Z') 



4r 



AS K {Z') 



where 



Here 



CY Lh (M Z >) 6/21 

0.985 7/6 



cJ^(m) = i + 



1 + 1.371 ^^ (1 



-in 



a s 

4-7T 



-2 log 



Mi 
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ft 



^ + - 
2 3 



4f 



(73) 



(74) 



(75) 



represents 0(a s ) QCD corrections to Z' tree- level exchange |53| and the two factors 
involving 



ocf\M zl ) 

(6)/ \ 

a s '{m t ) 



(76) 



represent together NLO QCD renormalization group evolution from mt to Mz' as given 
in |54| . The renormalization scheme dependence of this evolution is cancelled by the one 
of L . We now explain the origin of r. 

As the operator structure in 331 models is unchanged with respect to the SM, the 
B q and Bk parameters that have been factored out are the same as in the SM. The 
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same applies to the renormalization group evolution in the SM and 331 models between 
low-energy scales relevant for hadronic matrix elements and fit = 0(mt) entering the 



evaluation of 772 and tjb 49 that have also been factored out. Therefore the departure 
of f from unity is governed by the renormalization group effects between fit and fiz' = 
O(Mz'), absent in the SM, and the difference in matching conditions between full and 
effective theories, involving tree diagrams in the 331 models but box diagrams in the SM. 



These are represented by C^ LL (M^/) and the numerical factor 0.985 49 , respectively. 
The latter factor describes flavour universal QCD correction to So(xt) in the SM and is 
usually included in tjb and r\2 |49|. The coefficient 1.371, calculated in [54j, corresponds 
to the effective theory with six flavours / = 6. We assume tacitly that the new fermions 
S,D,T are heavier than Z'. Their inclusion into this formula would have a very small 
impact on f. Note that f is free from hadronic uncertainties. It is also flavour universal, 
as all the flavour dependence is already included in 772, ?/\b, B q and Bk that have been 
factored out. 
We find then 

f{M Z f = 1 TeV) = 0.985, f{M z < = 3 TeV) = 0.9534. (77) 



3.4 New Box Contributions 

In principle one would also have to include new box diagram contributions with new 
charged and neutral gauge bosons and new quark exchanges. Also charged Higgs par- 
ticles can contribute. However, without any calculation one can convince oneself that 
these contributions are negligible. Indeed, compared to the SM contributions the new 
box contributions are suppressed automatically by a factor M^/My < 0.006, where 
My > 1 TeV stands for a new gauge boson. In models with new LR operators, like left- 
right symmetric models this suppression could be compensated by enhanced hadronic 
matrix elements of new operators and QCD renormalization group effects. But in 331 
models only SM operator is present and such enhancements are absent. Another en- 
hancement could be present through enhanced values of the elements of the mixing 
matrix Vl- However, our analysis shows that in order to suppress sufficiently tree-level 
Z' contributions to AF = 2 processes, the hierarchical structure of Vl resembles the 
structure of the CKM matrix. 

Box diagrams are also suppressed with respect to Z' tree-level contributions by a 
loop factor like l/(167r 2 ) and two additional vertices. Therefore they cannot compete 
with the latter. 

In summary it is safe to keep only Z' contributions, which significantly simplifies 
the analysis as new box diagrams depend generally on new parameters, like masses of 
new gauge bosons and fermions and new mixing parameters, which are absent in Z' 
contributions. 



3.5 Basic formulae for AF = 2 observables 

Having the mixing amplitudes M| 2 at hand we can calculate all relevant AF = 2 ob- 
servables. To this end we collect below those formulae that we used in our numerical 
analysis. 

The Ki — Ks mass difference is given by 

AM K = 2 [S (M$ gM + 5tf (M^) Np ] , (78) 
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and the CP-violating parameter ek by 



c ^v&^ [ ^ M ^« + ^ M ^ ] ' (79) 



where 93, = (43.51 ± 0.05)° and n e = 0.94 ± 0.02 (25}[55| takes into account that ip e ^ f 
and includes long distance effects in Q^r^) and ^(M^). 

For the mass differences in the -B^ — B® systems we have 

AM, = 2 |(M« 2 ) SM + (Mj 2 ) Np | (q = d,s). (80) 
Let us then write [56] 

Mf 2 = (Mf 2 ) SM + (M? 2 ) Np = (M* 2 ) SM C Bq e 2 ^ , (81) 

where 

' e 2ip , (3 « 22° , (82) 



Aff 2 



SM 



M? 2 



SM 



(^! S 2 )sm = I(^i s 2 )smI e 2iA , Ps - -1° • (83) 
Here the phases f3 and j3 s are defined through 

Vtd = \Vtd\e- ip and y ts = -|F ts |er^ . (84) 

We find then 

AM, = (AM q ) SM C Bq , (85) 

and 

S^Ks = sin(2/3 + 2(p Bd ) , = sin(2|/3 s | - 2y B J , (86) 

with the latter two observables being the coefficients of sin(AM<jt) and sin(AM s t) in the 
time dependent asymmetries in B® — > ipKs and B® — > ificf), respectively. Note that new 
phases are directly related to the phases of the functions S q : 

2<p Bq = -0f . (87) 



At this stage a few comments on the assumptions leading to expressions in (86) are 
in order. These simple formulae follow only if there are no weak phases in the decay 
amplitudes for B^ — > ipKs and B® — > ip<f> as is the case in the SM and also in the LHT 
model, where due to T-parity there are no new contributions to decay amplitudes at tree 
level so that these amplitudes are dominated by SM contributions [45] . Similarly in the 
model discussed in the present paper there are no new contributions to decay amplitudes 
relevant for S^,k s and and the formulae given above apply. 



4 Effective Hamiltonians for AF = 1 Decays 
4.1 Preliminaries 

The goal of the present section is to give formulae for the effective Hamiltonians relevant 
for rare K and B decays that in addition to SM one-loop contributions include tree 
level contributions from the Z' gauge boson. While new penguin and box diagrams 
involving other gauge bosons and heavy quarks can also contribute, they are subleading 
with respect to Z' tree-level contributions and we will neglect them in what follows. 
Moreover, in contrast to Z' contributions they involve additional parameters and are 
more model dependent. 
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Figure 2: Tree level contribution of Z' to the s — > dvv effective Hamiltonian. 



4.2 Effective Hamiltonian for s — > dvv 

The effective Hamiltonian for s — > dvv transitions resulting from Z-penguin and box 
diagrams is given in the SM as follows 



ft 



\ K 

e //J SM 



2 

5SM 



E 

=e,/i,T 



where x ? - 



(K) 



*nni><0 + Ar } X(x t ) {s llx P L d){v nil P L vt) + /i.e. , (88) 



V* s V id and Vij are the elements of the CKM matrix. 



X^ NL (x c ) and X(xt) comprise internal charm and top quark contributions, respectively. 
They are known to high accuracy including QCD corrections 57-59] and electroweak 
For convenience we have introduced g| M that is defined in ( 72 ) . 



corrections 



60 



61 



In the 331 model (88) is modified by the tree-level diagram in Fig. [2] A straight- 

K 



forward calculation of the diagram in Fig. 2 results in a new contribution to 



TJVV 

H eff 



Combining this contribution with the SM contribution in (88) we find 



vv 
eff. 



K 



2 

5SM 



E 

=e,/i,T 



+ h.c. . 

where X(K ) is given as follow^ 



X{K) = X{x t )+ LK ' L 



Af(Z') Ai d (Z>) 



= \X{K)\e ie x. 



Here 



X{x t ) = T]X 



x t + 2 3x t - 6 



+ 



lnxt 



Vx 



0.994 



(89) 



(90) 



(91) 



_x t - 1 ' (x t - l) 2 

results within the SM from Z-penguin and box diagrams, r/x is QCD correction to these 



diagrams 57,62] when rat = nit(mt). 

The relevant couplings A^(Z') are defined in the Appendix [a] As all NP contributions 
have been collected in the term proportional to \\ K \ X^ NL (x c ) contains only the SM 
contributions as stated above. 



4 In this section in order to increase the transparency of basic formulae we use the notation Xg = X(K) 
and X q = X(B q ) and similarly for other functions. 
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Figure 3: Tree level contribution of Z' to the b — > s£ + £ effective Hamiltonian. 



4.3 Effective Hamiltonian for b 



dvv and b — > svv 

dvv and b 



In order to generalize the result just obtained to the case of b 
transitions only two steps have to be performed: 

1. All flavour indices have to be adjusted appropriately. 

2. The charm quark contribution can be safely neglected in B physics. 
The effective Hamiltonian for b — > qvv (q = d, s) is then given as follows: 

tvv 1 B. 



SVV 



effl 



9sm E KV tb X(B q )](qrP L b)(Pa»PL 



with 



+h.c. , 



X(B q ) = X(x t ) + 



Z') 



v t *v tb 



\X(B q 



(92) 



(93) 



.9smM 2 z ,_ 

Again all relevant A^ 9 entries in the 331 model can be found in Appendix [a| 

Note that the functions X(K) and X(B q ) presented above depend on the quark 
flavours involved, through the flavour indices in the A^(Z') = s, d, b) couplings and 
through the CKM elements that have been factored out. While in principle A^(Z') 
could be aligned with the corresponding CKM factors, this is generally not the case and 
the functions in question become complex quantities that are flavour dependent. This 
should be contrasted with the case of the SM and CMFV models where K, Bj and B s 
systems are governed by a flavour-universal loop function X{xt) and the only flavour 
dependence enters through the CKM factors. Consequently, as we will see below, certain 
SM-relations and more generally CMFV-relations will be violated in the 331 model. 



4.4 Effective Hamiltonian for b — > d£ + £ and b 



i£+t 



The effective Hamiltonian for b — > q£ + £ (q = d, s) transitions in the 331 model is a 
simple generalization of the SM one. For (q = d, s) we have 



n 



off 



i=A,V 



where 



Qa = (WuPLbXhV), Qv = {qiaP L b){hn) 



(94) 



(95) 
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and H e ff(b —> 57) stands for the effective Hamiitonian for the b — > 57 transition that 
involves the dipole operators. An explicit formula for the latter Hamiitonian will be 
presented in the next subsection. 

In the 331 model the second term in (94) receives contributions from the usual SM 
penguin and box diagrams and the tree-level Z' contribution shown in Fig. [3] that domi- 
nates NP contributions. Usually such tree-level gauge boson exchanges generate primed 
operators obtained from the ones in (95) by replacing Pj, with Pr. But in the 331 model 
similarly to the SM right-handed FCNC's are absent and primed operators are irrelevant. 
We find then for the Wilson coefficients 



C A 



where we have defined 



C v = Y(x t ) -4 sin 2 6 w Z(x t ) 
-Y(x t ) 



1 1 Af{Z')Af{Z') 



1 



9 2 SM M 2 Z , 
1 Af(Z')Af(Z') 



v t *v tb 



Slu M l 



v t *v tb 



(96) 
(97) 



Af[Z') = Af(Z') + Af(Z'), 



Af{Z') 



A^(Z') 



(98) 



Here Y(xt) and Z{xt) are SM one-loop functions, analogous to X(xt), that result from 
various penguin and box diagrams. The relevant entries have been collected in Appendix 
|A} In particular, 



Y(x t ) = tjy 



4 3x t log x t 



8 \xt-l (xt-1) 2 



Tjy = 1.012, 



where t]y is QCD correction to these diagrams |57[|62| when rat = mt{ m t) 
relevant entries have been collected in Appendix |A"} 
Introducing 



(99) 
All other 



Y{B q ) = Y{x t ) + 



Af(Z') 



v t *v tb 



\Y(B q 



ie q 
e y, 



Z(B q ) = Z(x t ) + 



1 



2A^(Z') 



Af(Z') 



v t *v tb 



4sin 2 ^ M*,g* M 
the Wilson coefficients Cy and Ca can be cast in a SM-like form: 

Cv = V t *V tb [Y(B q ) - Asm 2 9 w Z(B q )}, 

C A = -V%V a Y(B q ). 



Z(B„ 



(100) 

(101) 

(102) 
(103) 



The effective Hamiitonian for s — > d£ + l~ transition can be obtained directly from 
formulae given above by replacing q by K, appropriately changing the flavour indices 
and neglecting the contributions of primed operators. Then the functions Y(K) and 
Z{K) enter: 



Y(K) = Y(x t ) + 



Af{Z') 
M l9 2 su 



Ag(gQ 

v t * s v u 



\Y(K)\e* 



(104) 
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Figure 4: New magnetic penguin diagram contributing to B — )■ X s 7 wt/i intermediate Z' . q 
denotes a SM quark. 



Z{K) = Z(xt) + — l — 
4 sin a 



w 



2Aff(Z') 

M 2 z ,g 2 SM 



|Z(^)|e ? ^. (105) 



As seen in the Appendix, in the 331 model the coupling Ay 1 {Z') is strongly suppressed 
and by one order of magnitude smaller than A^(Z'). Consequently Z' contribution to 
the coefficient Cy is very small and NP enters b — > st + l~ and b — > d£ + £~ transitions 
dominantly through the coefficient C\, that governs the decays B S( i — > u + yT . In our 
numerical analysis we will discuss only these decays, but we have checked that NP effects 
in b — > s d£ + l~ transitions are within present experimental and theoretical uncertainties. 

In the case of the minimal 331 model, Ay l (Z') is by a factor of three larger than 
A^(Z') but as the NP contributions in this model are tiny, this change is irrelevant for 
all practical purposes. 

4.5 Effective Hamiltonian for the B — > X s ") Decay 
4.5.1 Preliminaries 



Very detailed analyses of B — > X s j decay in 331 models have been presented in [23 42 , 
where further references can be found. It has been found that the dominant NP contri- 
butions come from the Higgs sector, while the gauge boson contributions are subleading. 
As the Higgs contributions involve other set of new parameters than the ones involved 
in our analysis, there is no impact on our analysis from the constraints on Higgs con- 
tributions to B — > A s 7 [23]. Similarly the subleading contributions involving the new 
gauge bosons V°, V° and and the new heavy quarks involve new parameters and 
their contributions can be suppressed if necessary without any impact on our analysis. 

An exception is the Z' contribution in Fig. [4] which involves only SM quarks and the 
parameters that enter our analysis. While on the basis of [23] we expect this contribu- 
tion to be small, we present here a new analysis which, as far as QCD corrections are 
concerned, goes beyond the analysis of the latter paper. Even if finally this contribution 
will turn out to be negligible in the 331 models, the formulae presented below could be 
useful for other models. 

Adopting the overall normalisation of the SM effective Hamiltonian we have 

n cS (b -> s 7 ) = --^-V t * s V tb [C 7l (fi b )Q 7l + Csg^Qsg] , (106) 
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where fi b = 0(m b ). The dipole operators are defined as 
fn b s a (T^ u P R b a F^ v , Qsg = ^2 



= v^m b s a a^P R b a F^ , Q 8G = z^m b s a a^ P R T^bpG» v . (107) 



In writing (106) we have dropped the primed operators that are obtained from (107) by 



replacing P R by Pl- In the SM the primed operators (RL) are suppressed by m s /m b 



relative to the ones in (106). This is also the case of the 331 models. We have also 
suppressed current-current operators which are important for the QCD analysis. We 
will include these effects in the final formulae at the end of this subsection. 

The coefficients Ci(/x&) are calculated from their initial values at high energy scales 
by means of renormalisation group methods. We distinguish between SM quark contri- 
butions with the matching scale fit = 0(mt) and the Z' quark contributions with the 
matching scale fj,z> = O(Mz')- While in the LO approximation the results depend on 
the choice of the matching scale, the experience shows that taking as the matching scale 
the largest mass in the diagram appears to be a very good choice at LO. The choices 
made above follow this strategy. 

We decompose next the Wilson coefficients at the scale \x b = 0(m b ) as the sum of 
the SM contribution and the Z' contributions: 

CM = Cf M (fi b ) + AC?' ( M6 ). (108) 

We recall that for the SM coefficients at fxt = 0{mt) we have (xt = mf/M^) without 
QCD corrections 

«w - + 'l^-'iT' s (uo) 

4.5.2 Z' contribution without QCD Corrections 

A general analysis of neutral gauge boson contributions to B — )■ X s ^ decay has been 
presented in [63]. In addition to SM-like LL contribution from Z' we have a new LR one, 
where L (R) stands for the Pl (Pr) projector in the basic penguin diagram involving 
the s(6)-quark. In the 331 models the LR contributions to B — > X s ^ originate from 
flavour conserving bbZ' and ssZ' couplings which in addition to left-handed component 



have also a right-handed component as clearly seen in (63). 

In what follows we present the results for a contribution of a fermion / carrying 
electric charge —1/3 and having the mass rrif. As only SM quarks with irif Mz' 
contribute, we can set mf = 0, whenever it is justified. 

We first decompose the Wilson coefficients ACf at the fxz 1 scale as the sum of the 
SM-like LL contribution and a new LR one: 

ACf 7 W) = A LL C%M + A iR Cf 7 W) , 

(111) 

AC&Oiz') = A LL Qb(M + A LR Cl G (vz>) ■ 
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Adapting the general formulae of |63| to our notation and denoting by / the down- 
quark exchanged in the diagram we find 

\i>L r z> (n s 2 1 Ml ^ Al s *(Z>)A{ b (Z>) _ 2 Ml 4^33 
77 { ^' } ~ 9 g* Ml 2f - V * Vtb - ~ 27 Mj, V* V tb 

A LL Ci G (fi z ,) = -3A LL C*(» Z ,) , 

where the summation is over SM down-quarks. As this contribution is independent of 
quark masses the summation over / can be performed implying a very simple result, 
in particular when the d-quark contribution suppressed by additional mixing angles is 
neglected. 

For LR Wilson coefficients we find in the case of /3 = 1/ y/3: 

A LR r z> (l v _ 2 1 Ml v m f A{ S *(Z') AgOgQ ^ 2 s\ v* 32 v 33 

77 ^ Z ' ] 3 Ml 2-, mb V * Vtb ~ 27 M f, 1 - V/ s Vfc ' 

A iK Qb(/xz') = -3AX(MZ') , 

(113) 

where again the summation is over SM down-quarks. As LR contributions arise from 
the chirality flip on the internal line, in this case explict dependence on mj is present. 
However in accordance with other approximations we should keep only the b-quark 
contribution, which gives the final result. As 

C^(x t ) = -0.193, Cl^(x t ) = -0.096, (114) 

it is evident that without inclusion of QCD corrections, the Z' contribution to B — > 

is totally negligible for the values of M z > considered. We will now demonstrate that this 

is also the case after the inclusion of QCD effects. 

4.5.3 Final Results including QCD corrections 

In order to complete the analysis of B — > X s ~f we have to include QCD corrections which 
play a very important role in this decay. In the SM these corrections are known at the 



NNLO level 64 . In the LR model a complete LO analysis has been done by Cho and 



Misiak |65| and after proper modification we can use their results in our model . In this 



context the recent analyses 63 , 66 turned out to be very useful. 
We find then 

AC^(fi b ) = K 7 (fi z ,) ACf 7 W) + « 8 0i Z ,) ACioM + Ar nt (w) • (H5) 



The last contribution in ( |115 ) results from the mixing of new neutral current-current 



operators generated from the Z' exchange that mix with the dipole operators. The 
renormalization group analysis of this contribution is very involved but fortunately the 



LO result is known from 63 . Therefore adapting the formulae (4.16), (4.17) and (5.6) 



of this paper to our notation we find 



A C Z ™V) = kI la ^ LA C f 2 M + Y&A ^ LA C^z>), (H6) 



A=L,R A=L,R 
f=u,c,t,d,s,b 
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where 



and 



A LA C^ 2 



2 M^At(Z')A{ f (Z' 



9 2 M% 



v: s v th 



2 M^Af*(Z')A™{Z') 



9 2 Ml 



v t *v tb 



(117) 
(118) 



The diagonal couplings A^(Z') introduce additional parameters. 

Finally, k's are the NP magic numbers listed in Tab. [2] that is based on |63j which 
used a s (M z = 91.1876 GeV) = 0.118. They have been obtained for [i h = 2.5 GeV as 
used in the SM calculations. 

Using these formulae we find for M Z i = 1 TeV 



ACf 7 '(/x 6 ) 



0(1(T 5 ) 



(119) 



which is negligible when compared with the SM value of —0.353. Therefore we will not 
consider the B — > X s ^f further. 





1 TeV 


5 TeV 


10 TeV 


K 7 
Kg 


0.457 
0.125 


0.408 
0.129 


0.390 
0.130 


u,c 
K LL 

K LL 

K d 
LL 
s,b 

K LL 

H d 


0.057 
-0.003 
-0.057 
0.090 
0.147 


0.076 
-0.002 
-0.072 
0.090 
0.163 


0.084 
-0.001 
-0.079 
0.090 
0.168 


u,c 
K LR 

K LR 
^d 
^LR 
s,b 
K LR 
is-d 
^LR 


0.128 
0.012 
-0.025 
-0.092 
0.665 


0.173 
0.023 
-0.036 
-0.106 
0.865 


0.193 
0.028 
-0.041 
-0.111 
0.953 



Table 2: The NP magic numbers relevant for QCD calculations lblft 
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5.1 K + — > -k + vv> and K L 



7T UU 



Having at hand the effective Hamiltonian for s — > dvv transitions derived in Section 4.2 
it is now straightforward to obtain explicit expressions for the branching ratios B{K + — > 
ir + vi>) and B{Ki — > tt^vv). Reviews of these two decays can be found in 67-69). 
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The branching ratios for the two K — > ixvv modes that follow from the Hamiltonian 
in Section 4.2 can be written generally as 



7T 1 VV) 



Bv(K L 



f !mX eS \ 2 f RcX ca „ fin \ 



where 70 



(5.36 ±0.026) • 10 



-li 



V J ' 



k l = (2.31 ±0.01) • 10" 



-10 



and [58-60,70,71 



P C {X) = 0.42 ± 0.03. 
The short distance contributions are described simply by 

^cff = V t * s V td X(K) = V t * s V td \X(K)\ e w *. 



(120) 
(121) 

(122) 
(123) 

(124) 



where X(K) is given in (90). 

These formulae are in fact very general and apply to all extensions of the SM. The 
correlation between the two branching ratios depends generally on two variables 
and 6^ and measuring these branching ratios one day will allow to determine them 
and compare them with model expectations. The numerical analysis of both decays is 
presented in Section M 



5.2 B -> {X s ,K,K*}vv> 



Following the general analysis of [72| , the branching ratios of the B — > {X s , K, K*}vv 
in the 331 model can be simply obtained by removing the right-handed current contri- 
butions: 



where 



B(B -> Kvv) = B(B ->• Kuu)sm x q 2 , 
B(B K*vv) = B(B -»• K*uu)su x Q 2 
B(B -> X s uv) = B(B -»• X s vt?) sm x Q 2 

= \X{B S )\ 
6 X(x t ) 



(125) 
(126) 
(127) 

(128) 



with X(B q ) defined in (93) 



Evidently, the ratios of these branching ratios are equal to the corresponding ratios 
in the SM and all three branching ratios depend on NP only through the variable g. It 
should also be emphasized that the average of the K* longitudinal polarization fraction 
Fl, also used in the studies of B — > K*£ + £~, takes in the absence of right-handed 
currents a fixed value 

(F L ) = 0.54 . (129) 
Therefore, the measurement of Fl constitutes an important test of 331 models. 



We should remark that the expressions in Eqs. (125)-(127), as well as the SM re- 



sults in (130), refer only to the short-distance contributions to these decays. The latter 
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are obtained from the corresponding total rates subtracting the reducible long-distance 
effects pointed out in [73| . 



The updated predictions for the SM branching ratios are 72 - 74 
B(B -)■ Kvv) SM = (3.64 ± 0.47) x 10~ 6 , 
B(B -> K*vv)sm = (7.2 ± 1.1) x 10~ 6 , 
B(B -> X a vu)sM. = (2-7 ± 0.2) x 10~ 5 , 



(130) 



to be compared with the experimental bounds 75-77] 

B(B -> Kvv) < 1.4 x 10~ 5 
B(B -»■ iTVz/) < 8.0 x 10~ 5 



B(S -»■ X s i/z7) < 6.4 x 10" 



(131) 



5.3 B 



d,s 



We will next consider the two super stars of the LHCb, the decays B^ s — > that 
suffer from helicity suppression in the SM. This suppression cannot be removed through 
the tree level exchange of Z' boson. Assuming that the CKM parameters have been 
determined independently of NP and are universal we find 



B(B q 



B(B q 



Y(B n 



Y(x t ) 



(132) 



where Y(B q ) is given in (100). The numerical results are given in Section [8j 

The branching ratios B(B q — > fi + are only sensitive to the absolute value of 
Y(B q ). However, as pointed out recently in 30,31] in the flavour precision era these 



decay could allow to get also some information on the phase of Y(B S ) and we want to 
investigate whether in the 331 model this effect is significant. 

First as stressed in |30[|3l[ [78| when comparing the theoretical branching ratio with 
experimental data quoted by LHCb, ATLAS and CMS, a correction factor has to be 
included which takes care of AT S effects that influence the extraction of this branching 
ratio from the data: 



Here 



with 



B{B S -> fi + fi ) th = r(y s ) B{B S -»• fi + fi ) exp , r(0) = 1. 
r{y s ) 



1 ' /: 



Vs = r Bs 



1 + A\ v y. 



0.088 ±0.014. 



(133) 
(134) 
(135) 



The quantity -4^ r is discussed below. 

As stressed in 12] it is a matter of choice whether the factor r(y s ) is included in 
the experimental branching ratio or in the theoretical calculation, provided r(y s ) is not 



3 We follow here presentation and notations of 30 31 
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significantly affected by NP. Once it is measured, its inclusion in the experimental value, 
as advocated in [30], should be favoured as it would have no impact on the theoretical 
calculations of branching ratios that do not depend on Ar s . As in the SM and CMFV 



*4^ r = 1 31 and the factor r(y s ) is universal, it is also a good idea to include this factor 



in experimental branching ratio. In this manner various CMFV relations remain intact. 



If this is done the experimental upper bound in (174) is reduced by 9% implying a more 
stringent upper bound of 3.8 x 10 -9 . For the latest discussions see 30 31 and 79 80 . 

If a given model predicts *4^r significantly different from unity and the dependence 
of r(y s ) on model parameters is large one may include this factor in the theoretical 
branching ratio: 

B(B S fi^fi ) corr = — r (136) 

r{y 3 ) 

with the result for the SM [30j[3l] 

B(B S ->• M + /0^r = ( 3 - 5 ± °- 2 ) • 10_9 - ( 137 ) 

It is this branching that should be compared in such a case with the results of LHCb, 
ATLAS and CMS. However, as we will see below in the 331 models considered by us A\ T 
is very close to unity and it is more convenient to include this effect in the experimental 
branching ratio. 

What is interesting is that in addition to A\y it is possible to define a CP-asymmetry 
St+p- |3lj analogous to S^k s aricl S^. Both observables can be measured one day 
providing additional tests of NP models. 

The authors of |3TJ,|soj provide general expressions for A\ T an d as functions 

of Wilson coefficients involved. Using these formulae we find in the 331 models very 
simple formulae that reflect the fact that Z' and not scalar operators dominate NP 
contributions: 

A X AT = cos(20*' - 2<p B .), S£ +/i _ = sin(20*- - 2<p Bs ) (138) 

Both A^ r and are theoretically clean observables. 

In the SM and CMFV models 

^Ar = l, S£+„-=0, r(y s ) = 0.912 ±0.014 (139) 

independently of NP parameters. As we will see in the next section, the first and conse- 
quently the third quantity above is valid rather accurately in the models considered here. 
But for 331 the CP-asymmetry S s + _ will definitely be different from zero. Moreover, 
it will be of interest to investigate how is correlated with S^Ks anci ^W- 

While Ar^ is very small and yd can be set to zero, in the case of — > [i + one 



can still consider the CP asymmetry 5^+^- 80 , for which we simply find 



S* +IX - =sm(2e*<-2<p Bd ) (140) 

Even if this asymmetry, similarly to S s + _, is experimentally challenging, it will be 

M /* 

interesting to find its size in the 331 model. 
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5.4 K, 



The discussion of the NP contributions to this decay is analogous to B^^ — > ^ + H~ ■ 
Again only the SM operator (V — A) ® (V — A) contributes and the real function Y(xt) 
is replaced by the complex function in Eq. (104). Only the so-called short distance (SD) 
part to a dispersive contribution to Ki — > fi + n~ can be reliably calculated. We have 
then following (8l] (A = 0.226) 



B{K L m + /Osd = 2.08 • 10- 9 [P c (Y K ) + A 2 R t \Y{K)\ cos^ 



iH 2 



where Rt is given in (178) and 

\Vcb\ 



A 



3 K 
1 Y i 



Pc {Yk) 



1 



Pc(Yi 



K 



(141) 



(142) 



with P c (Yk) = 0.113 ± 0.017 [82]. Here (3 and j3 a are the phases of V t d and Vj s defined 
in @. 

The extraction of the short distance part from the data is subject to considerable 
uncertainties. The most recent estimate gives [83] 

B(K L -> m+m~)sd < 2.5 ■ 10" 9 , (143) 

to be compared with (0.8±0.1) • 10~ 9 in the SM 82 . The numerical results are discussed 
in Section [H 



5.5 B + 



T + V 



5.5.1 Standard Model Results 

We now look at the tree-level decay B + — > t + u which in the SM is mediated by the 
exchange with the resulting branching ratio given by 



B(B- 



-t^Jsm 



G 2 F m B+ m 2 T 



mz 



m, 



l P g+\V u b\ 2 T B + 



(144) 



Evidently this result is subject to significant parametric uncertainties induced in (144) 
by F B + and V u b- However, it is expected that these uncertainties will be eliminated in 
this decade and a precise prediction will be possible. Anticipating this we will present 
the results for fixed values of these parameters. 

In the literature in order to find the SM prediction for this branching ratio one 
eliminates these uncertainties by using AM<j, AM^/AM, and S^k s 



84 



85 



and taking 

experimental values for these three quantities. This method has a weak point as the 
experimental values of AM^ S used in this strategy may not be the one corresponding to 
the true value of the SM. However, proceeding in this manner one finds |85| 



B(B + -> t+^)sm = (0.80 ± 0.12) x 10" 



(145) 



with a similar result obtained by the UTfit collaboration 84 and CKM- fitters |86|. 



Now, the experimental world avarage based on older results by BaBar 87 and Belle 
was (89l 



B{B + -> tV 



exp 



(1.65 ±0.34) x 10" 



(146) 
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which is roughly by a factor of 2 higher than the SM value. Very recently, new results 
have been provided by BaBar Collaboration |90|: 



-)• T + U) 



(1.79 ±0.48) x 10" 



BaBar 



1-4 



B(B 

and by Belle Collaboration [9l] : 

B (B + -> r + i/)exp = (0.72 ±° f 5 ±° f 1 ) x 10 
implying a new world average provided by the UTfit collaboration 12 

(0.99 ±0.25) x 10~ 4 



Belle 



B(B + ->• t + u) 



cxp 



WA, 



(147) 



(148) 



(149) 



which is consistent with the SM. As we will see below the central value in (149) for 
central lattice input corresponds to \V U \,\ = 0.0040. 

The full clarification of the situation will be provided by the data from Super-B 
machines at KEK and Rome. In the meantime hopefully improved values for F B + from 
lattice and \V u b\ from tree level decays will allow us to make a precise prediction for this 
decay without using the experimental value for AM d . 

5.5.2 Result in the 331 Model 

In the 331 model there are no new treedevel contributions to B + — > t + u t . This means 
that the model favours a value of \V u b\ corresponding to the inclusive determination of 
this CKM element. We will return to this point in our numerical analysis. 



6 Breakdown of CMFV Relation in the 331 Model 



The presence of new sources of flavour and CP violation in the 331 model modifies the 
usual CMFV relations [IT] . Therefore it is useful to generalize the most interesting 
among these relations to include possible breakdown of them as follows [2j: 



AM d m Br , B d F% 



AM, 



m Bl 



Vtd 

Vts 



r(AM) 



(150) 



B(B d 



B(B -)■ X d vu) 
B(B X a vv) 



Vtd 



Vt 



ts 



r\vv) 



_ T(B d ) m Bd F B d 

B(B s ^v+fi-) r{B s )m Bs F 2 Bs 



Vtd 



V 



ts 



r(/x + /i ). 



(151) 



(152) 



The quantities r(AM), r{yv} and r([i + /j, ) are all equal unity in models with CMFV. 
In the 331 model they can be entirely expressed in terms of the master functions Q: 



r(AM) 



Sd 



r{uu) 



X(B d 
X(B S 



Y(B d ) 
Y(B S ) 



(153) 
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Eliminating (V^/V^l from the three relations above allows to obtain three relations 
between observables that are universal within the CMFV models. In particular from 



(150) and (152) one finds 92 



B(B S 
B(B d 



) 



B d t(B 8 ) 
B s r{B d ) 



AM S 



r(AM) 
r(/i+//~) 



(154) 



that does not involve Fb q and consequently contains smaller hadronic uncertainties than 
the formulae considered above. It involves only measurable quantities except for the ratio 
B s / B d that is now known already from lattice calculations with respectable precision 
1.05 ± 0.07. A review of other CMFV relations is given in 



17 



[93J|94| 

We also recall that within CMFV models there is a correlation between AM S 
and AM d as all these quantities are described by the same universal function So(xt). As 
a result, the enhancement of one of these observables implies uniquely the enhancements 
of other two. As we will see below also this correlation is broken in the 331 models. 



7 Correlations within 331 Model 

7.1 Classification 

The CMFV correlations between various observables are violated in the 331 model, in 
particular the correlations between the K, B d and B s systems are less stringent. Yet we 
would like to emphasize that the fact that NP effects originate dominantly from tree-level 
Z' exchanges with only left-handed quark couplings involved and leptonic Z 1 couplings 
are fully described by SM parameters implies still very stringent correlations between 
departures from SM and CMFV in three classes of observables. In particular there are 
interesting correlations between AF = 2 and AF = 1 observables in each class which 
imply important tests for this NP scenario. 

Indeed in the three meson systems, K, B d and B s , NP effects are governed by 

Af{Z'), Af{Z>) A h £{Z>), (155) 

respectively and by the value of the Z 1 mass which is obviously common to these three 
systems. As seen in Appendix |a| the three A^(Z') depend only on four new parameters: 

«13) S23, <*1, $2 (156) 

with S13 and S23 being positive definite and Si in the range [0, 2tt]. Therefore for fixed 
Mz'i the Z 1 contributions to all processes analyzed by us depend on only these param- 
eters implying very strong correlations between NP effects to various observables. The 
relations discussed below make it clear why the correlations in question are present. 

As we will see in the next section §13 and §23 are small numbers implying that £13 
and C23 are close to unity. Therefore there is a very clear structure in the dependence 



of various observables on the parameters in (156). In view of this structure we expect 
stringent direct correlations between observables in the following three classes: 
Class A: 



which depend only on §23 and S2 



AM S , B s ^fi + fi-, Ai r , B^X s uv, (157) 
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Class B: 

AM d , S^ Ks , B^fr, (158) 

which depend only on §13 and 6%. 
Class C: 



e K , K + ^it + vi>, K l ^tt°uu, K L ^H + H~ : K l ^tt £ + £~, e'/e, (159) 

which depend on S13, §23 and 62 — 61. 

While the processes and related observables in classes A and B are at first sight 
uncorrelated from each other, the fact that the processes in class C depend on all four 
parameters implies indirect correlation between all classes when experimental constraints 
are taken into account. 

7.2 Master Relations 

In view of a very simple structure of Z' contributions, it is possible to derive the relations 
between the modifications of the master functions Si, Xi and Yi for arbitrary @ 7^ v3- 
We will first list these general expressions and apply them to the case j3 = l/v3. Sub- 
sequently we will give them for /5 = \/3- Finally we will discuss the implications for 
phenomenology performed in the next section, where we will see the implied correla- 
tions between various observables in explicit terms. In what follows we will denote the 
modifications of master functions F{ coming from Z' contributions simply by AFi. 
We first find the relation between the Z' effects in AF = 2 master functions 



AS K _ M 2 z ,gl M 



AS d AS* 4f 



Sdl ryl\ I 2 n oa f HJr \ 2 



Af(Z 



A b £{Z')A b £*{Z') 



3.68 / M z > \ m 02 



3TeVy 



i-(i + ^y w ). (i6o) 



Here and in following equations we set \Vtb\ = 1 and C13 = C23 = 1 if necessary, f is given 



in (j77|). 

Next relations between the Z' contributions to AF = 1 functions Xi and Si functions 
are given by 



a q - ~ L "7 r = ~a q ^ ( "-iP- ) , )= ' -"'"H , (161) 



-.2 



AX q _ n Af(Z') _ ^ 0.085 {3TeV\ [l - (l + \/3/3)s 

w 



where 



Moreover, 



a d = l, a s = -l. (162) 

(163) 



AX K _ AX S 

JASk y/KS*' 

where with our definitions there is no complex conjugate on the l.h.s of this equation. 
Finally, there is the relation between Z' contributions to Xi and Yf 

AY = [1 - (1 - V3f})s 2 w ] 
AXi 1 _ (1 + V3f3)s 2 w ' 
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Applying these formulae for the case (3 = 1/ v3 we find 

AS K 2.55 / M z < \ 2 
AS d AS* s ~ ~T V 3 TeV J ' 

AX q 0.055 /3TeV 



AYi 



(165) 
(166) 



AJQ - L86 ' (167) 
AYj = AZi . (168) 



For the y/3 case we find 



AS K 0.27 / M z > \ 2 



AS d AS* " f V3TeVy 
AX q 0.023 /3TeV\ 



^ASJ q ^Tf \ M z , ) 
AYi 



(169) 
(170) 



AXi ™ 
AYi = s 2 w AZ t = AX, , i = K,d,s. (172) 

7.3 Implications 

While our statements to be made below are more general, in giving numerical examples 
we will concentrate on the 331 case. The master relations listed above have a number of 
implications which one can see even before a detailed numerical analysis is performed: 



The relation (165) implies that after the experimental constraints on AM S ^ have 
been taken into account, the effects in ek must be smaller than in AM S ^ for values 
Mz' < 3 TeV. Indeed, NP contributions in AM S]l j are constrained to be at most of 
10% of the SM values and therefore the absolute values of AS q can be at most 0.25. 
As ri w 1, we find then ASk < 0.16, which implies a correction of at most 6%. 
However, as we will discuss later on for much larger M z > these effects will increase 
significantly. 



The relation (166) implies that NP effects in the processes with vv in the final state 
are rather small for M Z ' = 3 TeV. Requiring again that AS q < 0.25 and taking 
into account that X(xt) ~ 1.5, we find that for M Z i = 3 TeV only effects of at 
most 5% at the level of branching ratios are expected for b — > svi> transitions. For 
K + — > tt + uv and Ki — > ifivv these effects are even smaller. For M z > = 1 TeV the 
relevant mixing parameters S13 and §23 have to be decreased by roughly a factor of 
3 to satisfy AM S ^ constraint. The inspection of the dependence on mixing angles 
shows then that this decrease of the latter is compensated approximately by the 
decrease of M z > in the case of K + — > tt^vv and Kl — > 7r°z/z/. However, it is 
overcompensated in the case of b — > svv transitions, so that the modifications of 
the branching ratios in this case by 15% are possible. 
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• The relation (167) implies that NP effects in decays with [i [i~ in the final state 
can be larger than in the case of vv. This is not only because of the factor —1.9 
in this relation but also because Y(xt) is smaller than X(xt): Y{xt) ~ 1. Thus 
for Mz 1 = 3TeV effects of 10 — 15% at the level of the branching ratio are still 
allowed in B SjC [ — > and these effects are expected to be increased up to 30% 

for Mz> = 1 TeV. Interestingly the NP effects in B SyC i — > fVfjT and B —> X s vv 
are anti-correlated. The suppression of B s ^ — > fi + implies enhancement of B — > 
X s vv and vice versa. Similar comments apply to K + — > tt + vv and Kl — > fi + . 
Unfortunately, due to small NP effects in decays with vv in the final states, these 
anti-correlations will be difficult to test. 

Finally let us note that once the constraints on AS q are taken into account NP 
effects for (3 = \/3 both in ek and in all rare decays are so small that it will be difficult 
to distinguish this model from the SM on the basis of flavour violation in meson decays. 
Therefore, we will not consider this case further. 



8 Numerical Analysis 

8.1 Recent Data and Anomalies 

Our prime motivation for a new analysis of the 331 models are the most recent data on 



Sjhd, and B St d — > \i \i decays from LHCb that read (95 - 97 



S M = 0.002 ± 0.087, S^f = 0.035 ± 0.002, (173) 

B(B S = (3.2±J;|) x 10~ 9 , B(B S -> /i + /0 SM = (3-23 ± 0.27) x 10^ 9 , (174) 

B(B d -> < 9.4 x 10" 10 , B(B d ^ fi + fi-) SM = {1.07 ±0.10) x 10- 10 . (175) 

We have also shown the most recent direct SM predictions for B(B q — > [i + ^) with new 
lattice input [79] . Almost identical results are obtained by applying the CMFV relations 
between B{B q — > fi + and AM q within the SM j2 , 18 . In quoting these results we did 



not include the correction r(y s ) but it has to be taken into account either in the theory 
or experiment when the data improve. 

While the SM still survived another test, from the present perspective S^^ could still 
be found in the 2a range 

- 0.18 < S n < 0.18 (176) 

and finding it to be negative would be a clear signal of NP. Moreover finding it above 
0.1 would also be a signal of NP but not as pronounced as a negative value. 

Concerning B s — > its branching ratio can still be enhanced by a factor of 2 

and finding it (5 — 6) • 10~ 9 would be a clear signal of NP at work. Also finding it well 
below the SM value is still possible. Similarly we should note that the upper bound on 
B(Bd — > ^ + /i~) is still one order of magnitude above the SM value. It could turn out 
after all that it is B d — > /U + /U~ and not B s — > that will most clearly signal NP in 

these decays. 

Another motivation for our analysis is a number of anomalies seen from the point 
of view of the SM that could become more relevant as the data and the lattice input 
improve. In this context it should be emphasized that because of some visible ek — S^k s 



tension 24-29 within the SM the pattern of deviation from SM expectations depends on 
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whether ek or S^k s is used as a basic observable to fit the CKM parameters. As both 
observables can receive important contributions from NP, none of them is optimal for 
this goal. The solution to this problem will be solved one day by precise measurements 
of the CKM parameters with the help of tree-level decays. Unfortunately, the tension 
between the inclusive and exclusive determinations of \V u b\ |^]and the poor knowledge of 
the phase 7 from tree-level decays preclude this solution at present. 

In view of this, it was useful already for some time [2] to set 7 ~ 70° in the ballpark 
of tree- level determinations and consider two scenarios for \V u b\. 

• Exclusive (small) \V u b\ Scenario 1: \ek\ is smaller than its experimental deter- 
mination, while S^k s is close to the central experimental value. 

• Inclusive (large) \V u b\ Scenario 2: \er\ is consistent with its experimental de- 
termination, while S^Ks is significantly higher than its experimental value. 

Thus dependently which scenario is considered we will ask the 331 model to provide 
constructive NP contributions to \ek\ (Scenario 1) or destructive NP contributions to 
S^k s (Scenario 2) without spoiling the agreement with the data for S^k s (Scenario 1) 
and for \ek\ (Scenario 2). 

This strategy turned out to be useful in the analyses of other models pj and helped 
to identify the correlation S^k s ~ ^i><t> ~ Wub\ i n the context of models with [/(2) 3 flavour 
symmetry [l4j. Yet, one should emphasize the following difference between these two 
scenarios. In Scenario 1, the central value of \ek\ is visibly smaller than the very precise 
data but the still significant parametric uncertainty due to | V^t, | 4 dependence in \ek\ and 
a large uncertainty in the charm contribution found at the NNLO level in 52 does not 
make this anomaly as pronounced as this is the case of Scenario 2, where large \V u b\ 
implies definitely a value of S^k s that is by 2 — 3a above the data. 

While models with many new parameters can face successfully both scenarios remov- 
ing the deviations from the data for certain range of their parameters it is not obvious 
from the beginning whether with improved data and new lattice input the 331 model will 
be successful in this respect, in particular when we will request from the 331 model to 
remove simultaneously other anomalies that we will list below. In fact in simpler models 
without new sources of flavour violation only one scenario for \V u b\ can be admitted as 
only in that scenario a given model has a chance to fit ek and S^Ks simultaneously. For 
instance as reviewed in [2] models with CMFV select Scenario 1, while the 2HDM with 
MFV and flavour blind phases, 2HDM MFV , selects Scenario 2 for \V u b\. 

Now the tensions within the SM discussed above constitute only a subset of visible 
deviations of its predictions from the data. In Table[3]we illustrate the SM predictions for 
some of these observables in both scenarios setting 7 = 68° . We only want to emphasize 
two points related to this table. 

First, the SM branching ratio for B + — > t + u t in Scenario 1 appears to be visibly below 
the data, although the latter are not very precise. Moreover, as discussed previously the 



most recent Belle result 91 in (148) and the world average in (149) are fully consistent 
with the SM. In Scenario 2 there is good agreement with the data. Until the experimental 
number will be clarified from the present perspective in Scenario 1 for \V u b\ models 
providing an enhancement of this branching ratio should be favoured. 

What is also striking in this table is that with the new lattice input in Table [4] the 
predicted central values of AM S and AM^, although slightly above the data, are both 



6 For a recent review see 98 . 
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Scenario 1 " 


Scenario 2: 


Experiment 


\ £ k\ 


1.72(26) • 10~ 3 


2.15(32) • 10~ 3 


2.228(11) • 10" 3 


B(B+ -> t + u t ) 


0.62(14) • 10~ 4 


1.02(20) • 10~ 4 


0.99(25) • 10- 4 


(sin 2/3) true 


0.623(25) 


0.770(23) 


0.679(20) 


AM, [ps- 1 ] 


19.0(21) 


19.0(21) 


17.77(12) 


AM, [ps- 1 ] 


0.56(6) 


0.56(6) 


0.507(4) 



Table 3: SM prediction for various observables for \ V u b\ = 3.1 ■ 10 3 and \ V u b\ = 4.0 • 10 3 and 
7 = 68° compared to experiment. 



in a good agreement with the latter when hadronic uncertainties are taken into account. 
In particular the central value of the ratio AM S / AMj is very close to the data. These 
results depend strongly on the lattice input and in the case of AM, on the value of 7. 
Therefore to get a better insight both lattice input and the tree level determination of 
7 have to improve. From the present perspective, models providing 10% suppression of 
both AM S and AM, with respect to SM values appear to be slightly favoured. As pointed 
out in [99] this is not possible within the models with CMFV and as demonstrated in |2j 
these models do not offer a good simultaneous discription of AM S( j and ek- It is then 
interesting to see whether the 331 model is performing better in this respect. 

8.2 Strategy 

As already advertised at the beginning of our paper, it is not our goal to present a full- 
fledged numerical analysis of the 331 model including present theoretical, parametric 
and experimental uncertainties as this would only wash out various correlations between 
various observables that we would like to emphasize. Recent sophisticated Monte-Carlo 
analyses in the context of other models and model independent studies pf|[8| [l0|[l00] 
that in particular used some kind of an average between inclusive and exclusive values 
of \V u b\i while certainly interesting, do not allow often for a transparent presentation of 
such correlations that are seen in our simplified approach. 

In view of the comments above our strategy in our numerical analysis should begin 
with finding out whether it is the inclusive or exclusive value of \V U \\ which is favoured 
by the 331 model. As our discussion of previous section indicates, the imposition of 
the experimental constraints on AM S( j for 1 TeV < Mz> < 3 TeV implies very small 
NP effects in ek- Consequently the 331 model favours a value of \V u b\ that is closer to 
inclusive determinations. Then not only ek is in good agreement with the data, but also 
the branching ratio for B + —> t + u t is closer to the experimental value than in the SM. 
Therefore we will perform our numerical analysis within Scenario 2 setting the CKM 
parameters at the following values 
Scenario 2: 

\V US \ = 0.2252, \V cb \ = 0.0406, 7 = 68°, \V ub \ = 4.0 • 10~ 3 . (177) 

The values for \V US \ and \Vd,\ correspond to their central values measured in tree level 
decays and the value of \V u b\ is in the ballpark of inclusive determinations. We checked 
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G F = 1.16637(1) x 10 
M w = 80.385(15) GeV 
sin 2 9 W = 0.23116(13) 
a(M z ) = 1/127.9 
a s {M z ) = 0.1184(7) 



-5 



GeV" 



m u (2GeV) = (2.1 ±0.1) Me V 
m d (2 GeV) = (4.73 ± 0. 12) MeV 

m s (2 GeV) = (93.4 ± 1.1) MeV 



(1.279 ±0.013) GeV 
4.19" 



m c (m c ) 

m b (m b ) = 4.191°;^ GeV 
m t {m t ) = 163(1) GeV 
M t = 172.9 ±0.6 ±0.9 GeV 



m K = 497.614(24) MeV 

F K = 156.1(11) MeV 

B K = 0.764(10) 

K e = 0.94(2) 

T] X = 1.87(76) 

r/2 = 0.5765(65) 

r] 3 = 0.496(47) 

AM K = 0.5292(9) x 10~ 2 

\e K \ = 2.228(11) x 10~ 3 



ps" 



89 
89 
89 
89 
89 



94 
94 

94 



101 



89 



94 , 102 



25 



B{B -> X sl ) 
B(B + -> r 



v) 



(3.55 ±0.24 ±0.09) x 
= (0.99 ±0.25) x 10~ 4 



10" 



t b ± = (1641 ±8) x 10" a ps 



89 



89 
94 
94 
55 
52 
49 
51 
89 
89 



89 
12 



89 



m Bd = 5279.5(3) MeV 
m Bs = 5366.3(6) MeV 
F Bd = (190.6 ±4.7) MeV 
F Bs = (227.6 ± 5.0) MeV 
B Bd = 1.26(11) 
B Bs = 1.33(6) 
B B jB Bd = 1.05(7) 



F Bn 



226(15) MeV 



F Ba \JB B3 = 279(15) MeV 

f = 1.237(32) 

Vb = 0.55(1) 

AM d = 0.507(4) ps- 1 

AM S = 17.77(12) ps^ 1 

S i , Ks = 0.679(20) 

S u = 0.002 ± 0.087 

t{B s ) = 1.471(25) ps 
r(B d ) = 1.518(7) ps 



: 0.2252(9) 
(40.6 ± 1.3) x 10" 



49 



89 
89 
94 
94 
94 
94 
94 



94 



94 
94 
50 
89 
89 
89 
95 



103| 
1031 



|K S | 
\V cb \ 

\V^ cL \ = (4.27 ±0.38) x 10- 3 

I "t/excl 
I V ub 



|89| 
|89| 
19\ 



(3.38 ±0.36) x 10- 3 g 



Table 4: Values of the experimental and theoretical quantities used as input parameters. 



that varying \ V u b\ within ±5% has only minor impact on our results. Other inputs are 
collected in Table HI 

Having fixed the four parameters of the CKM matrix, the "true" values of the angle 
j3 and of the element \Vtd\ are obtained from the unitarity of the CKM matrix: 

1 — Rb cos 7 
Rb sin 7 



\Vtd\ = \V U s\\V cb \Rt, Rt = Jl + R 2 b -2R b co S ~f, cot/3 = - j— - 1 , (178) 



where 

We find then 

\V U \ = 0.0399 , |Vfai| = 8.51 • 10" 3 (S2) . (180) 
Concerning the direct lower bound on Mz> from collider experiments, the most strin- 



gent bounds are provided by CMS experiment 104 . The precise value depends on the 
model considered. While for the so-called sequential Z' the lower bound for Mz> is in 
the ballpark of 2.5 TeV, in other models values as low as 1 TeV are still possible. In order 
to be on the safe side we will choose as our minimal value Mz> = 1 TeV and will also 
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provide values for Mz> = 3TeV. With the help of the formulae in subsection [7] it should 
be possible to estimate approximately, how our results would change for other values 
of Mz>- In this context let us notice that in view of moderate NP effects found in 331 
model the modifications of rare decay branching ratios due to NP will be governed by the 
interference of SM and NP contributions and consequently will be inversely proportional 
to Mz> as the corrections to the master functions. 



8.3 Simplified Analysis 

We first perform a simplified analysis of AM^ S , S$k s an d >SW m order to identify oases 



in the space of four parameters in (156) for which these four observables are consistent 
with experiment. To this end we set all other input parameters at their central values 
but in order to take partially hadronic and experimental uncertainties into account we 
require the 331 model to reproduce the data for AM S ^ within ±5% and the data on 
S-i/iKs an d S^cf, within experimental 2a. 

Specifically, our search is governed by the following allowed ranges: 



16.9/ps < AM, < 18.7/ps, -0.18 < < 0.18, (181) 

0.48/ps < AM d < 0.53/ps, 0.64 < S 4 , Ks < 0.72. (182) 

The search for these oases is simplified by the fact that the pair (AM S , S^) depends 
only on (S23, 62), while the pair (AM^, S^,k s ) on ly on (^13, <^i)- The result of this search 
for Mz> = ITeV is shown on the left in Figs. [5] and [6] for (§23, ^2) and (si3,5i), respec- 
tively. The red regions correspond to the allowed ranges for AM^ S , while the blue ones to 
the corresponding ranges for S^Ks anci S^. The overlap between red and blue regions 
identifies the oases we were looking for. Analogous plots for Mz> = 3 TeV are shown on 
the right in the same Figs. 5]and[6j From these plots we extract several oases that are 
collected in Tables [5] and [6 We denote by Ai(Mz') and Bi(Mz'), the oases for B s and 
Bd system, respectively. We observe the following pattern: 

• The increase of Mz< by a factor of three, allows to increases S13 and S23 by the same 
factor. This structure is evident from the formulae for AS q . However, in view of 



the relation (166) this change will have impact on rare B S; d decays, making the NP 



effects in these decays for Mz> = 3 TeV smaller. The inspection of the formulae 



for AXk shows that when the constraints from B-physics in (181) and (182) are 
imposed, the increase of Mz> is compensated in rare K decays by the increase of 
S13 and S23 so that these decays practically do not depend on Mz>- In the case of 
ek the same phenomenon even increases the room for NP effects with increasing 
Mz> for values of several TeV, where all FCNC constraints can be satisfied. 

• Simultaneously the ranges for Si remain unchanged. 

• For each oasis with a given Si there is another oasis with Si shifted by 180°. 

• The oases with i = 2, 4 are very small. 

In the rest of this subsection we will confine our numerical analysis to these oases, 
investigating whether some of them can be excluded by other constraints. In this context 
one should note that the oases in the (523,^2) and (si3,<5i) spaces are without other 
constraints independent of each other. Consequently when the full space of four free 
variables is considered we deal really with 16 oases. Fortunately, as we will demonstrate 
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Figure 5: Ranges for AM S (red region) and (blue region) satisfying the bounds in eq. 
(181). The plot on the left is obtained for Mz> = 1 TeV, that on the right for Mz> =3 TeV. 
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Figure 6: Ranges for AMd (red region) and S^k s (blue region) satisfying the bounds in 



Eq. (182). The plot on the left is obtained for Mz> = 1 TeV, that on the right for Mz> = 3 
TeV. 



at the end of this section, 8 oases involving at least one of the oases can be 

eliminated by £k and AM S /AM^ constraints alone. Among the 8 remaining ones, we 
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Table 5: Oases in the space (§23,52) f or Mz> = 1 TeV and Mz> = 3TeV. The sign of S s + 
chooses the oasis uniquely. The same applies to the pair and AB(B S — > as discussed 

in the text. 



will concentrate at first on the following four oases: 

(Ai,Bi), (A U B 3 ), (A 3 ,B 1 ), (A 3 ,B 3 ). (183) 

We will see that when all observables discussed by us will be measured we will be able 
to identify uniquely the oasis that has the best chance to describe the data correctly. 

The case of the small oases involving one of the oases B2 , B4 will be discussed in 
Section 18.51 

As a final comment, we observe that the oases reported in Tables [5] and [6] actually 
describe squares in the spaces (£23, £2) and (§13, Si), while the corresponding regions 
in Figs. [5] and [6] have more complicated shapes. Indeed, in our numerical analysis of 
the various observables we have varied the parameters in the true oases, requiring that 
constraints ( |181[ ) and (182) are satisfied. 



8.4 The Search for the Optimal Oasis 

We will now demonstrate that it is possible to identify uniquely one of the oases listed in 



( 183 ) by considering simultaneously a subset of particular observables considered in our 
analysis. Once the optimal oasis has been selected, the correlations between different 
observables in this oasis will tell us whether the 331 model is consistent with experiment. 

Before entering the search, it should be emphasized that in view of the small number 
of free parameters involved, the patterns of flavour violation and CP-violation depends 
in the crucial manner on two conditions which led to the oases in Tables [5] and [6j 

• The requirement that S^k s has to be suppressed significantly below its SM value 
to agree with data, 

• Both AM S and AM^ must also be suppressed below their SM values. 

While the first requirement is certain in view of small contributions to Ek and implied 
large value of \ V u b\, the second depends crucially on the non-perturbative parameters and 
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AB(B d -> fj+fj,-) 




Bi(l) 


0.0020 - 0.0032 
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£2(1) 


0.0078 - 0.0082 
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£3(1) 


0.0020 - 0.0032 
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+ 




B 4 (l) 


0.0078 - 0.0082 
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#i(3) 


0.0063 - 0.0099 
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£ 2 (3) 


0.024- 0.025 


92° - 95° 
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£ 3 (3) 


0.0063 - 0.0099 


308° - 330° 


+ 




S 4 (3) 


0.024- 0.025 


272° - 275° 







Table 6: Oases in the space (§13, 5i) for Mz> = 1 TeV and Mz> = 3 TeV. The enhancement or 
suppression of B(B d — > with respect to the SM value chooses the oasis uniquely. The 

same applies to the sign of ■ 



it will be of interest to monitor how the patterns of flavour violation presented below 
may change when the lattice calculations will improve in precision. 

Inspecting the expressions for various observables in different oases, we have identified 
the fastest route to the optimal oasis. We describe this route in the first step below. In 
the following steps we will demonstrate how the correlations between various observables 
can give additional tests once the analysis is confined to a particular oasis. 

Step 1: 

It turns out that the optimal oasis can be found by just measuring two observables: 

B(B d ^ M +/0, (184) 

that can be considered as two coordinates. Already the sign of shifts of them with 
respect to the SM values identifies the optimal oasis This we summarize schematically 
in Fig. [7J 

We observe: 

• (Ai,B\) oasis (blue) is chosen when 

s l+»- < °> B ( Bd -> < B{ - Bd -> M + /Osm- (185) 

• (Ai,Bs) oasis (green) is chosen when 

S* +M - < 0, B(B d -> (1+ /T) > B(B d -> m + M")sm. (186) 

• (As,B\) oasis (yellow) is chosen when 

S* +A( _ > 0, B(B d -> M +/0 < B(B d -> /i + /x-) SM . (187) 

• (A3, S3) oasis (purple) is chosen when 

S '^- > °> B ( B * -> ^"V~) > B ( Bd -> ^ + /Osm- (188) 



7 Also the sign of distinguishes between big oases. But this asymmetry is much harder to measure 



than B(Bd —> fx + u ). 
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Figure 7: B(B d — > fi + fi ) versus S* + . The plot on the left is obtained for M z < = 1 TeV, that 
on the right for Mz> = 3 TeV. The colours indicate the oases to which the parameters belong: 
(A 1 ,B 1 ): blue, (A U B 3 ): green, {A^B X ): yellow, (A 3 ,B 3 ): purple, (A 1 ,B 2 ): red, {A U B A ): 
gray, (A^,B2): magenta, [A^^B^): brown. Red point: central value of SM prediction. 




Figure 8: versus for Mz> = 1 TeV (left) and Mz> = 3 TeV (right). Oasis A\: blue 

(lower region), A 3 : purple (upper region). Red point: central value of SM prediction. 



In Fig. [7] we also display the case in which (S13, 5\) belong to one of the two small oases 
B2, -B4, but we discuss this case below in Section [83} 

In Fig. |8j we show S^ +f ^_ vs S^. Again the requirement of suppression of AM S 
requires Sf +I _ to be non-zero. A negative value of Sf +I _ chooses scenario A± (blue), 

MM MM 

while a positive one scenario ^3 (purple). Note that in both scenarios the sign of is 
not fixed yet but it will be fixed by invoking B{B S — > /U + /x~) below. What is particularly 
remarkable in this plot that for Mz' = 1 TeV, \S^ + | can reach values as high as 0.5 



when \S, 



ip(f> 



0.2. Smaller values are found for M Z ' = 3 TeV. 



In Fig. 9] we show S^k s vs B(Bd — > fi + fi Q The requirement on S^k s and AMd in 



(182) forces B(B<i — > /x + jj, ) to differ from the SM value but the sign of this departure 



depends on the oasis considered. Here distinction is made between B\ (yellow) and B3 



8 We should remark that the central values for B(B d [i + fi-) SM = 1.0 x 10 -10 and B(B S /i+^-) SM = 
3.1 x 10 -9 shown in the plots correspond to fixed CKM parameters chosen by us and differ from the ones 



listed in (174) and (175) but are fully consistent with them. 
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Figure 9: S^k s versus B(B d — > ji + ji ) for M z > = 1 TeV (left) and M z > = 3 TeV (right). Oasis 
B\. yellow (left region), B s : green (right region). Red point: central value of SM prediction. 



(green) for which B{Bd — > u + u~) is suppressed and enhanced with respect to the SM, 
respectively. These effects increase with decreasing S^k s - These enhancements and 
suppressions amount to at most ±25% and ±10% for M z > = ITeV and M z > = 3 TeV, 
respectively. 

The fact that S s + _ and B{B<i — > fi + fi~) are very powerful in identifying the optimal 
oasis can be understood as follows. S s + _ is governed by the phase of the function 
Y(B S ) that originates in the Z' contribution. It can distinguish between A\ and A3 
oasis because the new phase 62 in these two oases differs by 180° and consequently sin ($2 
relevant for this asymmetry differs by sign in these two oases. Calculating the imaginary 
part of Y{B S ) in (100) and taking into account that it is and not A^ that enters 
Y{B S ) one can convince oneself about the definite sign of S^ +fi _ in A\ and A3 oases as 
stated above. 

As far as B(Bd — > ^ + ^) is concerned, it is correlated with S^k s that is already 
well determined. Therefore, the range of 5\ cannot be large. B(Bd — > u + can then 
distinguish between B\ and B3 oases because cos S± differs by sign in these two oases. We 
find then destructive interference of Z' contribution with the SM contribution in oasis 
B\ and constructive one in oasis B3 implying the results quoted above. 

It should be noted that on the basis of AF = 2 processes such a distinction between 
these oases cannot be made because the relevant amplitudes are governed by 25± and 
26 2 which differ by 360°. 

We end this step by briefly discussing S^ + ^_ which as seen in Table |6j and Fig. 
can also help by means of its sign to distinguish between different oases, although 



becomes difficult for M z > = 3 TeV. We note also that the sign of _ is always 
opposite to the sign of the shift in the corresponding branching ratio, which can easily 
be understood by inspecting the ranges of S\. Moreover, the predictions for S d + _ are 
rather precise. For instance at M Z i = 1 TeV in oases B\ and B3, the allowed ranges 
for 5^+^- amount to [0.27,0.44] and [—0.38,-0.24], respectively. While being aware 
of the fact that this CP asymmetry will be very hard to measure, this example shows 
that it can provide in principle a very stringent test of the model in question and other 
extensions of the SM. We show this in Fig. 11 where we have also shown the results for 
small oases that we will discuss at the end of this section. 
Step 2: 
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0) 



0.80 




Figure 10: versus S^k s for Mz> = 1 TeV (left) and Mz< = 3 TeV (right). Oasis B\. 



yellow (upper region), B s : green (lower region). Red point: central value of SM prediction. 




Figure 11: B(B d -> versus S* + ^ for M z > = I TeV (left) and M z > = 3 TeV (right). 

Oasis B\: yellow, B 3 : green, B 2 : magenta, B^: brown. Red point: central value of SM 
prediction. 



The reason why B(B S — > /x + /i~) cannot be presently as powerful as B{Bd —> u + fi~) in 
the search for oases is the significant experimental error on 5^,^ with which this branching 
ratio is correlated. However, inspecting this correlation in a given oasis constitutes an 



important test of the model. We show this in Fig. 12 While in the oasis A\ (blue) 
increases (decreases) uniquely with increasing (decreasing) B(B S — > /i + ^ _ ), in the 
oasis ^3 (purple), the increase of 5^,^ implies uniquely a decrease of B(B S — > /i + /x _ ). 
Therefore, while B{B S — > alone cannot uniquely determine the optimal oasis, it 

can do it in collaboration with S^. Therefore, finding both these observables above 
or below their SM expectations, would select the oasis A\, while finding one of them 
enhanced and the other suppressed (opposite sign in the case of S^^) would select A3 as 
the optimal oasis. We indicate this pattern in Table [5] 

If the favoured oasis will be found to differ from the one found in step 1, the 331 
model will be in trouble. In terms of observables this correlation between Step 1 and 
Step 2 could be as follows. Let us assume that B(B S — > /i + /i _ ) will be found below 
its SM value. Then the measurement of will uniquely tell us whether A\ or A^ is 
the optimal scenario and consequently as seen in Fig. [8] and Table [5] we will be able to 
predict the sign of 5^+^- . Moreover, in the case of sufficiently different from zero, 
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Figure 12: versus B(B S — > n + jjT) for Mz> 
A±: blue (down left to up right), A 3 : purple ( 
of SM prediction. 
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1 TeV (left) and M z > = 3 TeV (right). Oasis 
left to down right). Red point: central value 



we will be able to determine not only the sign but also the magnitude of This 



discussion shows that we have a triple correlation — S^^ — B(B S — > ) in the 



331 model: once the sign of 5"*+^- is known a unique correlation — B(B S — > /i + u~) 
is found. If in addition one of these three observables is precisely known the other two 
can be strongly constrained. 

Finally, we also note that B(B S — > n + u~) can for Mz> = ITeV and M Z i = 3 TeV 
deviate from the SM value by ±30% and ±10%, respectively. 

Step 3: 

In Fig. 13 we show B(B,i — > fJ> + /J.~) versus B(B S — > n + for > 0.05 (upper- 



panel) and < —0.05 (lower-panel). We observe four different areas corresponding to 
the oases in question but the positions of these oases are different in these two plots. 
Step 4: 



In Fig. 14 we show B(B — > X s uu) vs B(B S —> ). This correlation is valid in any 



oasis due to the general relation (167) present in the 331 model. As expected, NP effects 
are significantly larger in B(B S — > fi + fi~) than B{B — > X s vv). 
Step 5: 



We have calculated B(Kl — > ir vv) vs B(K + — > -k^vv) (see Fig. 15). It turns out 
that in the oases [A\, B\) and (^3, -B3) (blue) both branching ratios are suppressed with 
respect to the SM values but they are both enhanced in scenarios (A\, B3) and (A3, B\) 
(green) . Combined with the remaining steps the measurements of these branching ratios 
could in principle contstitute an important test of the model in question. However, 
unfortunately, the deviations from SM expectations are at most 5% at the level of the 
branching ratios so that these correlations cannot be tested in a foreseeable future. 
On the other hand, finding experimentally both branching ratios significantly different 
from SM expectations would put the 331 model in trouble. As expected the NP effects 
basically do not depend on the mass of Z' . Also NP contributions to B(Kl — > /i + /x _ ), 
even if by a factor of two larger than in K + — > ^vv, are well within long distance 
uncertainties involved in this branching ratio. 



In Fig. 16 we plot AM, and AM^ as functions of \ek\ in order to see how the 



correlation between these three observables is modified within the 331 model. These plots 
show that this model does not suffer from AM s ^-|e^:| tension that has been identified 
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Figure 13: B(B d versus B(B S n + yr) for M z > = 1 TeV (left) and M z > = 3 TeV 

(right). In the upper panel: > 0.05, in the lower: < —0.05. Oases: (Ai,B{): blue, 
(Ai,Bs): green, (As,Bi): yellow, (A3, S3): purple. Red point: central value of SM prediction. 
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Figure 14: B(B -> X 8 vv) versus B(B S -> /or M Z / = 1 TeV (left) and M z > = 3 TeV 

(right). The two oases A\ (blue) and A3 (purple) overlap in most parts. Red point: central 
value of SM prediction. 



in CMFV models [2|. Indeed, imposing, as we do, that AM S ^ satisfy the experimental 
values, ek is within a few % from the experimental central value. As expected the room 
for NP contributions in ek is larger for larger Z 1 mass. 
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Figure 15: £(K L ->- tt°z/z/) versus vr+z/z/) /or M Z / = 1 TeV (left) and M z , = 3 TeV 

(right). (A\,B\) and (A 3 ,B 3 ): blue, (Ai,B 3 ) and (A 3 ,Bi): green. Red point: central value of 
SM prediction. 
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Figure 16: AM, and AM d versus /or M z , = 1 TeV (left) and M z > = 3 TeV (right). All 
oases (A\,B\), (Ai,B 3 ), (A 3 ,Bi), (A 3 ,B 3 ) overlap. Red point: central value of SM prediction. 



In Fig. 17 we plot A\ r vs S^. Only for Mz> = 1 TeV and significantly different 
from zero, does -A^r differ significantly from unity. 
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Figure 17: A^ T versus for M z > = 1 TeV (left) and M z > = 3 TeV (right). Oasis Ax: blue, 
A 3 : purple. 



8.5 A Brief Look at the Small Oases 

We will now briefly discuss the small oases. It turns out that they correspond to NP 
contributions to Mf 2 and M[ 2 that are roughly by a factor of two larger than their SM 
values but carry opposite signs. As AM^ S involve the absolute values of the mixing 
amplitudes these oases cannot be eliminated on the basis of them. The situation is 
analogous to the sign of the Wilson coefficient CV 7 that cannot be fixed by B — > X s j 
decay. To this end other observables have to be invoked and this is what we will do now 
in the case at hand. 

We anticipated that among the 16 possible oases obtained when the two pairs of 
parameters (§23, ^2) and (§13, <5i) belong to the regions denoted by Aj and B{ (i = 1, . . . 4) 
in Tables [5] and [6j respectively, it possible to rule out the 8 ones obtained when the pair 
(S23, 82) belongs either to A2 or to A4. This is a consequence of the requirement that ek 
does not differ more than 5% with respect to the experimental central value reported in 
Table |4j an uncertainty larger than the experimental range reported in the same table. 
This request is violated when (S23, 82) vary in A2 or in A4, allowing us to half the number 
of possible oases provided the theoretical and parametric uncertainties in ek could be 
lowered down to 5%. This is clearly not the case at present, primarly due to |V^,| 4 
dependence present in ek but could become realistic in the second half of this decade. 

In the previous subsection we discussed the case of the big oases, i.e. that in which 
the parameters belong to A4 and Bi with i = 1,3. There are still other 4 possible regions 
obtained when the parameters (si3,<5i) vary in B2 or B4. The resulting combinations of 
oases are therefore: 



(Ai,B 2 ), (A U B 4 ), 



(A 3 ,B 2 ), 



(A 3 ,B 4 



(189) 



At present we cannot rule out these regions, but from Fig. [7] and 11 we can observe that 



correlation between B(Bd — > u + u ) and S s 'f _ can help distinguishing among these four 



as well as for the big ones. In particular: 
• In (Ai,B±) (gray) 



MA* 



^W- < 0, 



B(B d ^n+n-)<B{B d 



)SM, 



(190) 
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• In {A U B 2 ) (red) 

S s ^. < 0, B(B d -> n+fT) > B(B d -> / u+/i-) S M, Sj +/1 _ > 0. (191) 

• In (^3,^4) (brown) 

S s ^ > 0, -> < Z3(B d -> m + /Osm, < 0. (192) 

• In (^3,^2) (magenta) 

S' +/1 _ > 0, BOBd -> /i+/0 > £(B d -> m"V)sm, ^ +At _ > 0. (193) 

However, compared to the big oases, these small oases are special since they give 
sharp predictions for the branching ratio of B d — > /x + fi~ decay and this may one day 
help to rule them out or confirm. We find: B(B d — > fi + fi~) = (1.74 ± 0.02) x 10~ 10 in 
(Ai,B 2 ) and (A 3 ,B 2 ) while B(B d -> = (0.911 ± 0.005) x 10" 10 in (A 1} B A ) and 

(^3,5 4 ), for M z > = 1 TeV, and B(B d -> M + A«~) = (1-172 ± 0.005) x lO" 10 in (Ai,5 2 ) 
and (A 3 ,S 2 ) while B(B d f*V) = (0.892 ± 0.001) x 10~ 10 in (Ai, J3 4 ) and (^ 3 ,-B 4 ), 
for M^' = 3 TeV. It can also be observed that, in the 3 TeV case there is no more overlap 
between the small and the big oases. 

Finally, we note that the asymmetry S d + _ is very precisely predicted in the small 

fj, 

oases. For instance at Mz> = 1 TeV in oases B 2 and B A , we find for S^ + ^_ = — 0.66±0.03 
and S* + = -0.38 ± 0.08, respectively. 



8.6 Deviations from CMFV Relations. 

The calculation of the various observables performed in the previous subsection within 
the 331 model allows us to discuss the departure from CMFV relations, i.e. the deviation 



of the functions in Eq. (153) from 1. We find that such functions vary in the following 



ranges when the parameters belong to the big oases: 
• 1 TeV 



r(AM) G [0.88,1.07] (194) 
r(vv) G [0.84,1.19] (195) 
r(fj, + fj,~) G [0.60,1.62] (196) 



• 3 TeV 

r(AM) G [0.88, 1.07] (197) 

r(vv) G [0.94, 1.06] (198) 

r(n + fT) G [0.85,1.18] . (199) 



We observe that deviation of r(AM) from 1 is approximately 10%, irrespective of the 
value of Mz>- r(yv) can deviate by almost 20% from 1 in the case Mz> = 1 TeV, while 
the deviation reduces to 6% for Mz 1 = 3 TeV. r(fi + fi~) is by far the function that mostly 
deviates from CMFV prediction: its values can differ from 1 by 60% for Mz' = 1 TeV, 
and almost 20% for M z > = 3 TeV. 
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Figure 18: Left: M z , = 1 TeV, right: M z > = 3 TeV. Oases: (A^B^: blue, (A U B 3 ): green, 
(A^Bx): yellow, (A 3 ,B 3 ): purple. 



The departure of these functions from 1 is anti-correlated in the case of r(u + u ) and 



r{yv\ as can be seen from Fig. 18 from which we argue that when r{yv) > 1 (< 1) 
it implies r(fi + u~) < 1 (> 1). The largest possible values of r{yu) are obtained when 
(si3,<5i) € Bi, while r(/i + /i~) is larger when (si3,£i) G B 3 . 

The case in which the parameters belong to the small oases is not displayed in Fig. 
18 However, we mention that in this case, for Mz> = 1 TeV, the values of r{yv) can be 
larger, while r(n + fi~) can reach even the value r(fi + fi~) = 2.2 in the region (A\,B2). 
For Mz> = 3 TeV the contribution of the small oases is almost indistinguishable from 
that of the big ones, except that r(fi + [i~) and r(v&) can be as large as 1.3 and 1.1, 
respectively. 

The violation of CMFV relations in the 331 model can be discussed also comparing 
the actual results for the branching ratios of B^^ —> to what one would obtain 

from Eq. (154) putting r = 1. This is done in Fig. 19 where we plot B{Bd — > fi + fj>~) 
vs B{B S —7- u + u~) and superimpose the band corresponding to the relation (154) with 
r = 1 and including the errors on the experimental determinations of AM S , AM^, 
t(B s ), T(Bd). On the other hand, we used B s /Bd = 1.05 ± 0.03, where the central value 
coincides with the lattice determination reported in Table |4j while we reduced the size 
of the uncertainty. An improvement in this direction would indeed substantially help 
understanding if and how much CMFV relations are violated in the 331 model. 



8.7 What if Z' is beyond the LHC Reach? 

Our numerical analysis assumed that Z' will be found at the LHC, that is Mz> is some- 
where in the range 1 — 3 TeV. In this range the B-physics constraints imply rather small 
values for §13 and S23 and consequently for their product that enters ek- NP effects in 
ek are then small and \V u b\ has to be sufficiently large in order that ek is consistent 
with the data. 

But what if Z 1 is beyond the LHC reach and its mass is as large as 10 TeV? We have 
analyzed this case with the following result: 



The increase of Mz> to 10 TeV still allows to satisfy the constraints in (181) and 
(182) provided the values of s\ 3 and §23 are appropriately increased (see below). 
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Figure 19: Left: M z > = 1 TeV, right: M z , = 3 TeV. Oases: (A^Bi): blue, (A-l,B 3 ): green, 
(AsjBi): yellow, (A^Bz): purple. Red point: central value of SM prediction. Pink band: 
CMFV relation from Eq. ( 154) with r = 1. 



• In turn, as already mentioned previously, NP effects in ek can be larger so that for 
Mz 1 = 10 TeV even the exclusive value for \V u b\ allows to reproduce the data for 
ek, which is not possible for the values of Mz> in the few TeV range. 

• NP effects in rare B-decays are strongly suppressed relatively to the effects pre- 
sented by us for 1 TeV < Mz 1 < 3 TeV so that here similarly to rare K decays 
everything is SM-like. 

Thus the main effect of NP in this case in the 331 model is to remove the anomalies 
in AF = 2 observables in the presence of exclusive value of \V u b\ without basically any 
impact on rare B and rare K decays. This rather uninteresting vision could become 
reality if eventually B(B + — > t + u t ) will turn out to be SM-like, favouring exclusive 
value of \Vub\. 

We have investigated this case in more detail by setting \V u b\ = 3.1 • 10~ 3 , Mz> = 
10 TeV and imposing additional constraint: 

0.75 < - - . < 1.25, 2.0 x 1(T 3 < \e k \ < 2.5 x 10" 3 . (200) 

This additional constraint is necessary as now NP effects in ek for S23 and S13 allowed 
by .B-physics AF = 2 observables in the presence of Mz> = 10 TeV can be large. We 
find then approximately: 

0.05 < s 2 3 < 0.12, 0.016 < S13 < 0.030 (201) 

and the ranges for <5i 5 2 slightly shifted to lower values. The reason for some departure 
from the expected M z > dependence in the range discussed until now is the additional 



constraint in ( 200 ) which is immaterial for Mz 1 of a few TeV. 

The fact that maximal values of ek increase with M Z i when only B-physics con- 
straints are taken into account is related to the increase of S23 and §13 as already ex- 
plained. But when the A sd (Z') reaches its maximal value, further increase of M z > will 
also decrease NP contributions to ek so that for very large Mz> SM results are obtained 
as required by decoupling of NP. However, this happens only for Mz 1 > 3000 TeV. This 
result is just a confirmation of the known fact that if the FCNC Z' couplings to quarks 
are 0(1), £k puts a very strong constraint on the scale of NP. 
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8.8 Sensitivity to Non- Pert ur bat ive Parameters 

In the present paper we have used the most recent lattice input to find out that the SM 
values for AM S and AM^ were visibly above the data, when the hadronic and CKM 
uncertainties were reduced down to ±5%. Here we would like to emphasize that some 
of the correlations presented by us would change if AM S and AMj would both be found 
below the data. In order to illustrate this point we have reduced the values of Bs q by 
20% so that now 

(AM s )gM = 15.2 ps- 1 , (AMrf)sM = 0.45 ps" 1 . (202) 

This result would be welcomed by CMFV models in which AM S and AM^ can only be 
enhanced keeping their ratio fixed. However, the enhanced value of S^k s in the large 
\V u i,\ scenario still would be problematic for these models. 

The 331 model faces this new situation as follows. It chooses the phases 8\ and 5 2 , 
dependently on the oasis considered close to 0° or 180°. For system we find: 

-12°<5 2 <12° (green), 167° < 62 < 192° (yellow) (203) 

with 0.0023 < S13 < 0.0036. For B s system we find: 

- 30° < 5 2 < 35° (purple), 150° < 6 2 < 215° (blue) (204) 

with 0.01 < S23 < 0.019. 

In both cases AM S and AM^ being sensitive to 25i are enhanced. In the case of 
B® — Bg system, where the asymmetry is small both in the SM and data, the 
apperance of a small 25 2 is rather natural. In the case of B^ — B® system the addition 
of a new contribution with a new phase much smaller than the SM phase lowers S^k s 
through the interference of both contributions, as desired. 

The question then arises what impact does this new structure have on correlations 



presented above. Setting Mz> = 1 TeV we answer this question in Figs. 20 and 21 for 
most interesting correlations. We observe: 



As seen in Fig. 20 in the B^ meson system, the general structure of correlations is 



unchanged but the NP effects this time in B^ — > are found to be much larger 



than in S + _ 



As seen in Fig. 21 the impact on the B s system is much larger. Now the roles 



of B s — > /i + jj, and 5*^+^- in the search for optimal oasis are interchanged. The 
suppressions or enhancements of B{B S — > relatively to the SM uniquely 

identify the optimal oasis. Moreover, B(B S — > is bound to be different from 

the SM values. 5^+^- can now be vanishing and can only identify the optimal oasis 
in combination with S^^, moreover it is now significantly smaller than previously. 

This analysis demonstrates clearly that it is crucial the have precise determinations 
of flavour observables within the SM in order to be able to search indirectly for NP. 



9 Summary 

In this paper we have reconsidered the flavour structure of the 331 models concentrating 
the phenomenology on a particular model, the 331 model. In this model the dominant 
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Figure 20: S$k s versus B{Bd — > H + f> ) (left) and B{Bd —> fi + fi ) versus (right). Red 

point: central value of SM prediction. For colour coding see (203). 
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Figure 21: S*+ versus (left) and versus B(B d — > ) (right). Red point: central 
value of SM prediction. For colour coding see (204). 



NP contributions come from Z' tree-level exchanges. For Mz' of order of few TeV, visible 
departures from the SM expectations are possible. As a detailed analysis of phenomeno- 
logical implications and related correlations have been presented in the previous section 
we only list few highlights: 

• Favouring the inclusive value of |V^|, for ITeV < Mz> < 3 TeV, the 331 model 
agrees well with the most recent B + — > t + v t data, removes the ek — S^k s tension 
present in the SM and the ek — AM St d tension present in CMFV models. Simul- 
taneously, while differing from the SM, it is consistent with the present data on 
B s ,d -> and S^. 

• We identify four big oases in the parameter space that can be uniquely distinguished 
by the values of B(Bd — >• and . Also the asymmetry is helpful 
in this respect. 

• In each oasis there is a definite correlation between B(B S — > and S^^ so 
that also these two variables, when considered simultaneously can give valuable 
information about the model. 
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• There are also four oases that cannot be ruled out at present but as predictions 
for various observables are very sharp in these oases, it is likely that they will not 
survive more precise data. 

• A characteristic feature of the model is the anti-correlation between B(B S — > fi + fj>~) 
and B(B X s vv) . 

• Z' contributions to b — > svv transitions, K + — > ir + vv and Kl — > ifivv are found 
typically below 5 — 10%. 

• The violation of CMFV relations are clearly visible with the largest one related to 
B s ,d -> decays. 

• We have investigated the case of Mz> = 10 TeV, demonstrating that for such high 
values of Z' mass the pattern of departures from the SM changes, allowing to 
remove all anomalies even for exclusive values of \ V u b\- But then NP effects in rare 
B and K decays are very small. 

We are aware of the fact that some of the correlations presented by us would be 
washed out if we included all existing uncertainties. Yet, our simplified numerical anal- 
ysis had as the main goal to illustrate how the decrease of theoretical, parametric and 
experimental uncertainties in the coming years might allow to exhibit certain features of 
NP, even if the deviations from the SM will be only moderate. In fact within the coming 
years the size of the assumed uncetainties in our analysis could likely become reality. 

The correlations between various observables found in this paper will allow to distin- 
guish it from other simple secenarios. In particular we have demonstrated that sizable 
violation of CMFV relations are still possible in this model. Here we would like to make 
a brief comparision of the 331 model with the 2HDM MF y and MU(2) 3 models. The 
most striking differences are as follows: 

and 331 model share the property that the NP contributions to 
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2HDM MW 

ek are small favouring thereby inclusive value of \V U \\. In 2HDM MFV , the data on 
S^k s imply then automatically a positive value of > 0.15 |2|. In 331 model this 
is not required and any value within the LHCb range is still possible. Concerning 
AF = 1 transitions let us just mention that in 2HDM MFV there is some tendency 
that with increasing S^a, also the lower bound on B(B S — > fj> + fi~) increases. In 



331 model, as seen in Fig. 12, the correlation — B(B S — > /j ', n~) is much more 
transparent. Moreover, while in the oasis A\ these two observables are correlated, 
they are anti-correlated in the oasis A3. In both models B{B S — > fj, + [i~) can be 
smaller and larger than the SM value, but in 2HDM MFV the deviations can be 
larger, a characteristic property of scalar currents. 

In the case of MU (2) 3 models, the most striking differences from the 331 model are 
the correlation S^k s ~ ^i>4> ~ \ Vub\ 14 that is absent in the 331 model and the fact 



that MU(2) 3 models can be made consistent with the FCNC data for exclusive value 
of \V u b\i while this is rather difficult in the 331 model for Mz' < 5 TeV. Finding 
in the future that exclusive value for \V u b\ is chosen by nature and Mz> < 5 TeV 
would favour MU(2) 3 models over 331 model. On the other hand if the inclusive 
value of \V u b\ is the right one and S^^ is found to be negative, this would put 
MU(2) 3 models into difficulties while the 331 model would face such data without 
any problems. 

In view of these definite findings we are looking forward to improved experimental 
data and improved lattice calculations. The plots presented by us should facilitate 
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monitoring the future confrontations of the 331 model with the data and to find out 
whether this simple model can satisfactorily describe the observables considered by us. 

The generalization of our analysis to Z 1 models with both right-handed and left- 
handed couplings and sizable NP contributions to ek and rare K decays is presented 
in 
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A Appendix 

We collect the formulae for the couplings of Z' to quarks and leptons especially for 
= v3 and ^ that enter 
Case 1: arbitrary f3 



/3 = \/3 and ^ that enter various expressions in our paper. 



9 



^l(Z') = , 



Af(Z' 



g [l-(l + V3/3)s 



2V3c w Jl-(l + f3 2 )s 



w 



Af(Z') = Af(Z') , 
Af(Z') = 



c w ^l-(l+P 2 ), w 



2~ ' 



Case 2: (3 = ^3 

Af{Z') 



A b L d (Z>) 



A b £{Z>) 



Af(Z') 



9 


V3cw\ 


/> 


- 4s 2 
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\/3cvka 


/> 


- 4s 2 
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i{5 2 -Si) 



Cw v h v 3i = 1-84 gsi 3 s 2 3Ci 3 e 



c w v h v 3i = -1.84 5Si 3 ci 3 C23e tSl 



V3c w Jl-4s 2 w 

0.090 g. 



-1.84 gs 2 3C 13 c 2 3e 
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Af (Z') = ^Af (Z') = Af(Z') = Af(Z') , (206e) 

Af(Z') = 0.271 g . (206f) 

Case 3: B = 4r 

v 3 

Af(Z') = 9 (? W vt 2 v zl = 0.61 <7Si 3 S23Ci 3 e i(<52 - 5l) , (207a) 

A b L d (Z') = 9 c 2 w v* 33 v 31 = -0.61 gSnhz^e-^ , (207b) 

A b £{Z') = 9 c 2 w v* 33 v 32 = -0.61 gs 2 3C 2 13 C23e-^ , (207c) 

Af (Z') = 9 ^"jf^ L = 0.213 5 , (207d) 

Af (Z') = Af (Z') , (207e) 

Af(Z') = 9 S y = 2S ^ 2 Af(Z') = -0.183 g (207f) 

Af (Z') = -0.396 5 , (207g) 

Af (Z 1 ) = 0.030 g. (207h) 
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